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CHAPTER 1 . INTRODUCTION 


During the last decade, the development of millimeter 
wave astronomy has dramatically improved our understanding 
of the nature of the interstellar medium. This has 
resulted primarily from observations of the spectral line 
radiation from rotational transitions of interstellar 
molecules. More than 50 molecular species have been 
identified in a few very bright sources whose physical 
conditions enable excition of the molecules [68]. Giant 

molecular clouds have been discovered in the Galaxy. 

) 

Molecular clouds have been found to constitute an 
appreciable fraction of the interstellar medium by mass 
and to be the raw material from which stars evolve. 

Studies of these molecular clouds and their distribution 
in space have led to an improved understanding of star 
formation processes [16], stellar mass loss [85] and 
galactic structure [9]« The relation between molecular 
clouds and supernova remnants [26], and between molecular 
clouds and young star clusters [56] are two among many 
current topics of investigation. 

Although a great deal of work has been done, many 
worthy experiments, such as large scale mapping of 
extragalactic sources, and searches for new interstellar 
molecules, become feasible only with instruments of 
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significantly greater sensitivity than those which 
presently exist. It is the main purpose of this research 
project to develop a new kind of receiver, using a 
superconductor - insulator - superconductor quasiparticle 
tunnel junction as the mixing element, to overcome the 
limitations inherent in all previously existing 
millimeter-wave receivers. 

As the radiation from interstellar molecules is 
extremely weak, sensitive spectral line receivers are 
needed to detect these very faint signals. Heterodyne 
receivers are commonly used in radio-telescopes to achieve 
high spectral resolution, which is necessary in order to 
determine the velocity of the molecular clouds. The 
heterodyne receiver consists of a "mixer” followed by an 
IP amplifier. The receiver noise temperature is: 


T 

1 R " 


-M 


+ LT 


IF 


( 1 . 1 ) 


where L is the mixer conversion loss, T M is the mixer 
noise temperature and Tjj, is the equivalent noise 
temperature of the IF amplifier. The two terms on the 
right side of eq. (1.1) are usually of comparable 
magnitude, so it is clear that to achieve a major 
improvement in T^, both mixer parameters, L and T^ must be 
improved . 
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For a heterodyne mixer receiver operating in the 
Rayleigh-Jeans ' limit, i.e. hv<kT, the sensitivity is 
given by the radiometric equation [69]: 


6T 


kt r 

( T Of ) 1 ' 2 


0 - 2 ) 


where «T is the smallest detectable signal, t is the 
observation time required to detect <5T, K is a constant 
and Af is the pre-detection bandwidth. From this equation 
it follows that the observating time required to detect a 
small signal is proportional to the square of the 
receiver's noise temperature. In order to exploit fully 
the potential of the telescope, it is necessary to improve 
the receiver performance to the point where the overall 
system noise temperature is comparable to or smaller than 
the sky noise (typically 135 K at 115*3 GHz). 

In the past decade, several different kinds of mixers 
have been developed, and GaAs Schottky barrier diode 
mixers have been the most commonly used for millimeter 
wavelength radio astromony. A cooled Schottky diode 
receiver, with about 125 K receiver noise temperature, 
operating at 90 - 120 GHz, has been reported [47]* 
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However, since this technology is very well developed, 
there is little prospect for any further dramatic 
improvement in the performance of semiconductor diode 
mixers. liquid helium cooled InSb hot electron bolometer 
mixers [45] have also been successfully operated, but 
their performance is severely limited by their very narrow 
instantaneous bandwidth (</■ 1 MHz). 

Several different kinds of superconductive tunneling 
devices have also been extensively studied for use in 
heterodyne receivers, among them are the point-contact 
Josephson junction, superconducting micro-bridge, 
super-Schottky diode and SIN junction. Point-contact and 
micro-bridge Josephson junction mixers belong to the same 
category, in which the extremely nonlinear Josephson 
tunneling currents are used to achieve mixing. Although 
these devices can achieve conversion gain, their intrinsic 
high noise properties [67] have eliminated them from the 
competition. Furthermore, due to their lack of mechanical 
stablity, experimental results with point-contact 
junctions are usually very hard to reproduce and the 
junction's thermal recyclability is also very poor. 

Unlike the Josephson effect mixer mentioned above, 
the super-Schottky (superconductor-semiconductor) diode 
mixer and the superconductor - insulator - normal metal 
(SIN) junction mixer, use the nonlinear quasiparticle 
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current-voltage characteristic to achieve mixing. The 
super Schottky diode mixer is limited by its series 
spreading resistance and junction capacitance, to 
relatively low (i.e. microwave) frequencies [60]. The 
performance of an SIN array mixer at 73-5 GHz has been 
reported [53] and the result is very encouraging. 

However, the I-V characteristic of the SIN junction is not 
as strongly nonlinear as that of the SIS device, and the 
resistance of the normal electrode may also tend to 
degrade the performance. SIN mixers are not, therefore, 
the best choice for the 100 GHz range. Since SIN mixers 
are intrinsically free from Josephson effect noise, they 
are very promising for use at higher frequencies (perhaps 
300 - 1000 GHz). 

Recently, the superconductor - insulator - 
superconductor (SIS) quasiparticle tunnel junction has 
been demonstrated by Richards and Shen [49] , Dolan, 
Phillips and Woody [14], and Rudner and Claeson [51 ] to be 
useable as a low noise millimeter wave mixer. It is known 
that there are two current carrying mechanisms in an SIS 
tunnel junction: the Josephson (pair) tunneling current 
and the quasiparticle (single electron) tunneling current. 
In an SIS mixer, heterodyne mixing is performed using the 
extremely nonlinear quasiparticle branch of the I-V 


characteristic which has a sudden onset of conduction at 
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the full superconducting gap voltage 2 a / e arising from the 
effective energy gap in the density states of the 
superconductors. From an ideal B-C-S calculation it is 
shown that this I-V characteristic approaches that of an 
ideal resistive diode near the gap voltage [22, 80]. 
Therefore, according to classical mixer theory, an ideal 
SIS quasiparticle mixer should have very good mixing 
performance (e.g. for an image rejecting mixer, L 0 dB) . 
However, both theoretical and experimental studies of the 
high frequency performance of this device have indicated 
that, when the photon energy, fto/e, at the input frequency 
becomes comparable to the voltage scale of the I-V curve 
nonlinearity, the classical theory is no longer able to 
describe the performance. Nonclassical phenomena such as 
conversion gain in a resistive mixer, and quantum-limited 
sensitivity, have been predicted by the "quantum mixer 
theory", recently developed by Tucker [ 73 ] and later 
verified by experiments. Recently, two other nonclassical 
phenomena, negative output impedance and infinite 
available gain, both predicted by Tucker's theory, have 
been observed in this device. Several complete SIS 
receivers, operating from 45 GHz to 250 GHz with 
performance comparable to Schottky diode receivers, have 
been reported by different groups [3>4, 40, 41 ,46,64,65]* 
With their potentially low noise properties and very low 
pump (10) power requirement, it is likely that the SIS 
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receiver will be the prime choice in the future for radio 
astronomy in the millimeter wave frequency range. 

This research is directed towards studying the 
performance of the SIS quasiparticle tunnel junction 
receiver in the 110-118 GHz region of the radio spectrum, 
with the ultimate objective of constructing an 
ultra-sensitive SIS receiver for 115*3 GHz CO line 
observations. Although Tucker's quantum mixer theory has 
provided the theoretical framework for analysis and design 
of the SIS mixer, this theory is rather complicated in its 
derivation and is very hard to interpret compared with the 
familiar classical theory. Furthermore, only under 
certain approximations (e.g. the 5-frequency 
approximation) is it possible to obtain a tractable set of 
equations which can be used to characterize mixer 
performance. Even these simplified equations contain 
infinite series of Bessel function products and a highly 
nonlinear junction response function. Hence, most 
theoretical calculations have hitherto been done by 
extensive computer simulation, and have been difficult to 
interpret. Apart from these dif f iculities , the practical 
situation is even more complicated. In this work, we have 
therefore taken a primarily experimental approach. In the 
following chapters we shall thoroughly discuss the design, 
fabrication, and optimization of our SIS quasiparticle 
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mixers. The performance of these mixers is in excellent 
agreement with the theoretical predictions using Tucker's 
theory assuming the 3-frequency approximation. This work 
has led to the construction of an SIS receiver which is 
now successfully operating on the Columbia-GISS CO Sky 
Survey telescope on the roof of the Pupin laboratory. It 
is the most sensitive receiver ever built in this 
frequency range, with a single sideband receiver noise 
temperature of 83 K at 115*3 GHz. Observations of the ^CO 
emission from molecular clouds beyond the solar circle at 
t = 80°, made using this SIS receiver, will be described. 

The success of this receiver has not only 
dramatically improved the observation efficiency of the 
Columbia-GISS telescope, but has also demonstrated that 
the SIS mixer has the potential to become the preferred 
front-end in future millimeter wave radio astronomy 
receivers . 

The following outline describes the contents of each 
chapter of this thesis, and indicates the original 
contributions : 

Chapter 2 contains a detailed description of Tucker's 
quantum mixer theory (including mixer noise theory), which 
is the fundation of this research project. We then apply 
this theory to the SIS mixer in two special cases: 1) 
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equal signal and image terminations (i.e. the DSB case), 
and 2) the image frequency short-circuited (SSB case). 

The conditions under which it is possible* to achieve 
negative output impedance and infinite available gain are 
investigated. A computer program was developed to predict 
and analyze the performance of an SIS mixer. 

Chapter 5 describes the design and fabrication of the 
junctions used in this work. The choice of the various 
fabrication parameters of SIS junctions for optimum mixer 
performace is discussed. 

In chapter 4, we describe the design of the mixers 
used in this work. Measurement techniques used to 
characterize the mixer mount are discussed. A microwave 
equivalent circuit of the mixer mount, including the 
effects of conductor loss, is derived, and used to study 
the available embedding impedances at the signal and image 
frequencies. This information is necessary for an 
accurate analysis of the performance of the SIS mixer, and 
for optimization of the stripline circuit in the next 
design iteration. 

Chapter 5 describes the experimental apparatus used 
for measuring the performance of the SIS mixer, and 
discusses the calibration of the measuring system. 


In chapter 6 we describe the mixing experiments in 
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which Tucker’s 3-frequency quantum mixer theory is 
quantitatively verified for the first time. In these 
experiments, the first observation of negative dc 
differential resistance and infinite available gain in an 
SIS mixer were made. These are quantum effects which are 
not possible in classical mixers. The best mixer 
performance we have achieved at 115-3 GHz are an SSB 
conversion loss of 6.9 dB and an SSB mixer noise 
temperature of 20 K, with a lower-sideband rejection > 25 
dB. 


The optimized SIS mixer described in chapter 6 is 
used to construct an ultra low-noise 110-118 GHz SIS 
receiver. Chapter 7 describes the construction and 
performance of this receiver, which is the most sensitive 
ever reported for this frequency range. 

In collaborating with Y.-L. Huang, we undertook some 
astronomical observations to demonstrate the power of this 
new receiver on the Columbia-GISS 4-ft telescope. The 
typical integration time per scan is only about 1 minute 
for an rms noise of 0.25 K with a bandwidth of 250 KHz per 
filterbank channel, demonstrating that the system is in 
fact 25 times faster than with the previous Schottky diode 
receiver. Observatioins of the distant molecular clouds 
towards the Cygnus-X region are presented in chapter 8. 
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In chapter 9 we discuss future planned work on SIS 
receivers and their ultimate limits. A summary of the 
present work is given. This concludes the thesis. 
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CHAPTER 2. THEORY 


2.1 Introduction 


Ever since its discovery, the superconducting 
(Josephson) tunnel junction has appeared likely to be 
useful as a very low noise amplifier, mixer or direct 
detector in the millimeter wavelength range. Its strong 
nonlinearity at low currents leads to extremely low shot 
noise, and it requires a very small amount of local 
oscillator power (at most a microwatt) to operate as a 
mixer. Also, the cryogenic operating environment suggests 
that a superconducting detector or mixer will have low 
thermal noise. In addition to the potential of these 
devices in fields such as millimeter wave astronomy, 
atmospheric physics and space communication, the 
development of the theory of the superconducting tunnel 
junction mixer and detector has its own academic interest, 
since this involves the study of the interaction process 
between the tunneling electrons and the input signal field 
(photon) . 

The basic principle of frequency conversion using 
nonlinear resistive elements has been extensively 
developed and interpreted by many authors [24,54,70,77] 
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employing a classical approach which assumes an 
instantaneous device response to the incident 
electromagnetic field. In classical theory, the small 
signal performance of a periodically pumped nonlinear 
device is usually described by a square admittance matrix. 
The elements of this matrix are simply related to the 
Fourier coefficients of the conductance and capacitance 
waveforms produced in the nonlinear device by the local 
oscillator. The small-signal current and voltage 
components at different sidebands will therefore be 
linearly related by this conversion matrix, and the device 
performance can be calculated in a straightforward way. 

When the input photon energy in volts, -ftu/e, is 
sufficiently large compared with the voltage scale of the 
device’s dc I-V nonlinearity, the classical theory is not 
applicable and the recently developed quantum mixer theory 
[73] must be used to calculate the performance of a mixer. 
This theory, which includes the quantized nature of the 
charge carriers in the nonlinear single-particle tunneling 
device, was originally constructed by Tucker 
[ 71 , 72 , 73 , 74 , 75 ] and was elaborated later by Feldman [ 17 ], 
Shen [ 57 ] and. Smith and Richards [ 62 ]. Tucker used 
photon-assisted tunneling theory to evaluate the elements 
of the mixer conversion matrix. Once this step has been 
done, the mixer analysis proceeds exactly as with the 
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classical theory. Noise properties of the tunnel junction 
mixer are also described by the quantum mixer theory. 

When this theory is applied to the SIS quasiparticle mixer 
with the input photon energy comparable to the voltage 
scale of junction's dc I-V nonlinearity, the predicted 
results differ strikingly from the classical analysis. 

For instance: reactive mixing terms appear in the mixer 
admittance matrix, and it becomes possible to achieve 
gain, even infinite available gain, under certain 
circumstances. Negative output impedance and quantum 
limited sensitivity are also possible. All of these 
(except possibly quantum limited sensitivity) have now 
been verified by experiments. Since this theory provides 
the theoretical framework for the design and analysis of 
our SIS quasiparticle mixer, it is outlined in some detail 
in this chapter. 

2.1.1 Basic Mixer Terminology 


In this section, we introduce some basic mixer 
terminology. The mixer contains one (or more) nonlinear 
elements such as a semiconductor Schottky barrier diode or 
a superconducting tunnel junction. The (weak) signal, at 
frequency u , and a (strong)' local oscillator wave, at 
frequency u , are coupled into these nonlinear devices, 
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where they generate, among other frequencies, an 
intermediate frequency (ujp) which is equal to the 
difference between signal and local oscillator 
frequencies. Other unwanted higher harmonics and 
sidebands are generally filtered out by the circuit. The 
performance of the simplest type of mixer is characterized 
by two quantities: the equivalent input noise temperature 
(Tjyj) and the mixer conversion loss (L g ) which is defined 
as : 


available power at the signal frequency 

L s = (2.1) 

power delivered to IP load 


Usually, a mixer will produce an IP output at u^p 

from either of two input frequencies : the signal, at u = 

s 

Up ± “jp, or the image, at - u^^ (The choice of + 

(jjp or 0 ^ - Ujp as the "signal" usually depends on the 
particular experimental circumstances). Clearly there are 
two separate conversion losses, L g and , at the signal 
and image frequencies. Under certain special conditions, 
equal powers may be applied to the mixer at both the 
signal and image frequencies, in which case it is useful 
to define a double-sideband conversion loss, L^, which is 
related to the two single-sideband quantities by [56]: 
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The quantity l c is appropriate, for example, when the 
input signal to the receiver is from a blackbody radiator, 
as is used for calibrating the receiver's noise 
temperature . 

A mixer for which L << L. (image rejecting) is 
called a single sideband (SSB) mixer, and its noise 
temperature, as measured using a blackbody radiator, will 
be denoted as T M OOTI . If I = L. , the mixer is called a 

ri , ooD S 1 

double sideband (DSB) mixer, and the noise temperature 
measured using a blackbody radiator is designated as 
Tjyj For some applications, a single sideband noise 

temperature is quoted for a DSB mixer. It is then given 

by [32,56]: 

T M,SSB = T M,DSB ^ + L s/ L i^ ( 2 - 3 ) 

Thus Tjyj ggg is twice T^ ^gg for a mixer which responds 
equally to both signal and image. It is known from the 
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classical theory of mixers [31,54] that the conversion 
loss of an ideal resistive mixer can be as low as 0 dB when 
there is complete image rejection, and 3 dB for a DSB 
mixer. Such an ideal mixer would generate no noise of its 
own, giving T M = 0 K. 

For spectroscopic observations, the spectral line 
appears only in one sideband, and it is more appropriate 
to use the single sideband values of noise temperatures 
and conversion loss. Thus in this thesis, except where 
explicitly noted, all numbers quoted are single sideband 
quantities . 
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2 . 2 Classical Mixer Theory - Generalized Heterodyne 
Receiver Model 


2.2.1 Introduction 


The fundamental concepts of classical mixer theory- 
have been studied in depth by Torrey and Whitmer [70], 
Uhlir [77], Saleh [54], and Held and Kerr [24]. Since 
this theory is well developed and is the essential 
framework for the quantum mixer theory, we shall summarize 
it in this section. 

In classical theory, the nonlinear mixer diode is 
represented by an equivalent circuit shown in Pig. 2-1 . A 
conversion admittance matrix can be formed to relate the 
small signal-sideband currents and voltages at the diode. 
The elements of this matrix are simply related to the 
Fourier components of the conductance and capacitance 
waveforms produced in the nonlinear diode by the local 
oscillator drive. The conversion loss and input and 
output impedances can then be determined from this 
admittance matrix provided that the embedding impedance of 
the mixer mount is known at various sideband frequencies. 
However, in performing the theoretical calculation of the 
mixer performance, it is necessary to know the correct 
waveform of the local oscillator voltage Vj^t) at the 
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diode and the resultant diode current. This problem can 
be accomplished by an appropriate large signal analysis 
[23,58]. For the moment we shall just assume that V^ 0 (t) 
is known and outline the small-signal mixer theory. 

2.2.2 Frequency Notation 

Saleh's frequency notation [54] will be used here. 

If a mixer is pumped at a frequency u^ and has an 
intermediate frequency u^, the frequencies resulting from 
a distortion-free mixer are those sideband frequencies 
which can be put in the following form: 

u n = nu p + u 0 * ^ = ~°° » **•* — 1 *0,1 ^ (2.4) I 

Thus u p and u^ represent the upper sideband, lower 
sideband, and sum frequencies respectively. Note that 
u n <0 for n = -®,...,-2,-1. The sideband frequency index n 
is also used as a subscript with various electrical 
quantities. For instance, Vp Vq, and V_^ represent the 
voltages at frequencies u^ , uq and u_^ ; Ip Iq and I_^ 
represent the currents at corresponding sidebands. 
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2.2.3 Conversion Admittance Matrix 


Assuming that a do voltage Vq, a local oscillator 
voltage V 10 (t) and a small signal v slg (t) are impressed 
across the nonlinear device, then the total voltage 
waveform is: 


V(t) = V 0 + 




W 10 


(2.5) 


Under the small signal approximation, the resultant 
current in the diode is given hy: 


<I(t)> = I LQ (t) 


i sig^ t ^ 


( 2 . 6 ) 


The small signal voltage and current can be expressed 
in terms of the individual sideband frequency component: 


V sig (t > - 


Re X 


m=— ® 


v e 
m 


iut 
m 


(2.7) 
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1 sig^ t ^ 


Re 2 


m= — 


1 rr, e 

m 


m 


( 2 . 8 ) 


Torrey and Vhitmer [70] have shown that for small 
signals these current and voltage components are linearly 
related by a matrix equation: 


m 




(2.9) 


the square matrix Y is called the conversion admittance 
matrix. For an ideal diode, the elements of this matrix 
are given by Torrey and Vhitmer and are: 


\m' 


= G_ 


m-m 1 


]u C , 
0 m m-m 1 


( 2 . 10 ) 


where G k and C k are the Fourier coefficients of the diode 
conductance G(t) and capacitance C(t): 


G(t) = 2 G k exp(jku t) , G k = G* k 

k=-® p 


( 2 . 11 ) 
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C(t) = 2 C k exp(jku t) 

k=-“ 



( 2 . 12 ) 


As shown in Fig. 2-1, the matrix Y is the admittance 
matrix of the intrinsic diode and can be regarded as a 
multifrequency multiport network, each port of which 
represents one sideband frequency, let and Y ffl be the 
source current and admittance termination at each port. 
Then an augmented network outlined by the broken line in 
Fig. 2-1 can be formed to represent the whole mixer which 
is characterized by the augmented admittance matrix Y' 
defined by: 


Y' 

mm 


f 


mm ' 


+ Y 6 , 

m mm ' 


and we have 


(2.13) 


I = i + Y v 
m m mm 


= S Y' ,v mt 


m 


, mm ' nr 


(2.14) 


inverting this matrix equation, we have 


v = y Z* , I , 
m mm m 

m 


(2.15) 


where 
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Iiz’ m ,n = iitI.ii " 1 

* + Vm.,'1" 1 (2 ' 16) 

In particular, the IP output voltage can he written 


as : 


0 


= 2 


m ' 




m ' 



( 2 . 17 ) 


where 


Om 


= Z 6m 


z oo 


( 2 . 18 ) 


and 


-eff 


= I 

m 


, *0m 


,I m' 


( 2 . 19 ) 
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can be regarded as an equivalent source current at the IP 
frequency connected in parallel with the load and the 
output admittance of the mixer (1 /Z^q). 

2.2.4 Mixer Conversion Loss 


The mixer conversion loss, from sideband u . to the IP 

uq, is defined as the ratio of the power available from 

the source Y. at the input port j to the power delivered 
J 

to the load Y^ at port 0. The available signal power at 
the port j is: 


P in 




( 2 . 20 ) 


Where I . is the source current at the input port j. The 
S J 

power delivered to the IP load Y^ is: 


P 


out 



0 


2 


( 2 . 21 ) 
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= ; g l'V 2|i sj 


where G. and G T are the real parts of Y. and Y T 
3 -h 3 L 

respectively. The mixer conversion loss is then given by: 


L 0 - 


P in 


out 


1 

4G j G i |z 6jl 2 


( 2 . 22 ) 


2.2.5 Mixer Port Impedances 


The impedance Z , which is defined as mixer's m-th 

port impedance (see Pig. 2-1), can be found by open 

circuiting the corresponding port termination (Y = 0) and 

measuring the impedance at port m of the augmented 

network. Therefore the port impedance is given by the 

mm-th element of the Z' matrix defined in (2.16) with Y = 

m 


0. That is: 
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Z mm,Y =0 
m 


I iy' , i rl n 

mm ' mm , 0 


(2.23) 


Where the subscript "0" indicates that the mm-th element 
of this matrix is evaluated with Y =0. For instance, 
the upper sideband input impedance of the mixer is: 


'1 , in 


1 1 

Y 10 

T 1-1 

01 

y oo +y l 

Y 0-1 

-11 

Y -1 0 

Y -1-1 +Y 


(2.24) 


'IF, out 


impedance 

• 

is : 

• 

• 

• 

• 

• 

.. Y n +Y 1 

Y 1 0 

Y 1-1 

Y 01 

O 

O 

Y 0-1 

•* Y -11 

Y -10 

Y -1-1 +Y 


(2.25) 
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2 . 3 Theory of Superconductive Tunneling 


2.3*1 Introduction 


A superconductor - insulator - superconductor (SIS) 

tunnel junction consists of two bulk superconducting 

electrodes separated by a thin insulator layer 

0 

approximately 20 ^ 40 A thick. Much has been learned 

about the nature of tunneling in the 20 years since 

Giaever's experiments on non-Josephson tunneling [20,21] 

and Josephson's theoretical treatment of superconducting 

tunnel structures that display phase coherence [29]. The 

goal of this section is to determine the tunnel junction 

quasiparticle current response in the presence of the 

applied field. Since we are only interested in 

non-Josephson tunneling, the Josephson effect (the 

supercurrent tunneling) will be neglected in this 
* 

treatment . 


2.3*2 Dynamic Tunneling Model 


* In experiments, this situation can be realized by 
applying a magnetic field to suppress the observable pair 
tunneling current, or by biasing the junction close to the 
gap voltage. 
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Following the treatment in Tucker [73] and Rogovin 
and Scalapino [50], the total system Hamiltonian for the 
tunnel junction in the presence of external fields is 
given by: 


H = H R + H 1 + eV ( t ) N i + H T 


(2.26) 


Where and are the full many-body Hamiltonians for 

the right- and left-hand electrodes. Coupling between the 

external field and the tunnel junction is characterized by 

the time-dependent potential eV(t) multiplied by the 

left-side number operator with = £ C k C k * 

k 

Hj is the transfer Hamiltonian which describes the 
coupling between the two electrodes and is to be 
considered as a perturbation. Following Cohen, Falicov, 
and Phillips [7], H^ can be written as: 



(2.27) 
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Where and are one-electron operators representing 
Bloch states on the left- and right-hand sides of the 
barrier. corresponds to the first term, in which an 

electron is transfered from a state q on the right to a 
state k on the left, while performs the reverse 
function. 

The tunnel matrix element T k + ( ^ + is given by: 


ktqt 


■s 2 f .. * 


q* * 


( 2 . 28 ) 


Where (* ) is the single particle wavefunction for a 

state in the left-hand (right-hand) metal assuming 
infinite barrier, ds is the area unit vector of the 
dielectric interface. The integral can be evaluated 
anywhere within the oxide. 

The junction current operator is given by: 


I 


«R 

= e 

dt 


le 


= — [h t , N r ] 
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16 4 . 

= — [Hi - HI] 

•k 

= I + + I_ (2.29) 

where I transfers a particle from right to left and I is 
its Rermitian conjugate. 

The tunnel junction is assumed to be in thermal 
equilibrium. The density matrix of the uncoupled system 
is determined by the junction Hamiltonian : 


e BH 

, Z = T e“ 6H (2.30) 

Z 


Where H = H^ + H^ - uN. The time-dependent potential 
term eV(t)N^ in the system Hamiltonian (2.26) is assumed 
only to modulate the energy of electrons on the left-hand 
side and to be turned on adiabatically . The addition of 
this term to the left side Hamiltonian is offset by the 
change in left side chemical potential: 


u L (t) = n R + eV(t) 


(2.31 ) 
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so that, the ensemble remains unaffected. In the 
Heisenberg picture, the single particle operator C^(t) on 
the left can be written as follows: 


t 

c k ( "t ) = cg(t) e xp[- ^ ^dt'V(f)] (2-32) 

Where C®(t) is the single particle operator in the 
absense of the applied field V(t). 

Separating the dc bias V Q from V(t), we can express 
the remaining time-varying contribution as Fourier 
components : 


ie 

exp{ 

1i 


t 

J dt'[V(t') - V 0 ]} 


r -i«'t 

/du'WU’Je (2.33) 


“00 
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then ( 2 . 32 ) becomes: 


03 

0 ieV n t f -iut 

c k (t) = Cj(t) exp[ -H-] /do'W( 0 ')e (2.34) 


The single particle operator on the right-side is just: 


c q (t) = c“(t) 


(2.35) 


The current operators (2.29) in the presence of the 
applied field can then be written as: 


1 (t) 


i°(t) 


exp[- 


leY O t ] 

ft 



-iu * t 
)e 


( 2 . 36 ) 


with 
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I°(t) 


- 

ft kq 


* r 0+ 
kq C q 


(t)C°(t) 


(2.37) 


and 


i + (t> = i+(t) 


( 2 . 38 ) 


From linear response theory [50], the current 
operator to the first order in the tunneling interaction 
is : 


I (1) (t) = I°(t) 



0 

/ 


dt 1 [h 2( t ' ) , I°(t)]e(t-f) 


0 


(2.39) 


Since for the uncoupled system <I U (t)> = 0, the 


expectation value for the total current can be written as 
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-U dt ' 


<I(t)> = - - / dt'<[l(t) , H T (t’ )]> 0 e(t-t' ) ( 2 . 40 ) 


CD 

2 f 

— Im I dt ' i 


ie(t-t')<[l + (t), I_(t’)]> ( 


Using (2.36), we get 


2 r * -i( u '- u ")t 

<I(t)> = - Im / d 0 'd 0 "W( 0 ' )W ( u " ) e 


co 

J dt'i e (t-t' )<[l°(t), I°(t)]> Q 

— co 


x exp[-i( u'+eV 0 /-h) (t-t * ) ] (2.41) 


Define the retarded current commutator evaluated in the 
equilibrium system as: 


X + _(t-f) = ie(t-t')<[l°(t), I°(f)]> 0 


( 2 . 42 ) 
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Then its Fourier transform is given by: 


(u) =— I d(t-t')X. (t-t' )exp[-iu(t-f )] (2.43) 

2tt J + " 


Let us define the complex response function j(u) as: 


4 IT 

d (u) = — X + _(u) (2.44) 


Then the average current (2.41) can be written as: 


a> 

/ * -i((j , -u")'t 

du • du"¥( u ' )¥ ( u ")e 

— CO 


x j ( u ' + eV^/fi) 


(2.45) 


If the junction is only biased by a dc voltage (i.e. V(t) 
= Vq), then the Fourier components of the phase factor 
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(2.33) are given by: 

W(u') = 6 (u • ) (2.46) 

then equation (2.45) becomes: 

<I(t)> = wv 

= Im[j(eV 0 /h)] (2.47) 

Where I^ c (Vq) is junction's dc current-voltage 
characteristic. 

The real part of j(u) is related to the imaginary 
part of j(o) by the Kramers-Kronig relation [81 ] : 


* This is a direct consequence of the causality 
principle. In the classical treatment, the current 
carriers crossing the junction barrier are assumed to 
respond instantaneously to the applied external field. 
However, at higher frequencies, i.e. approaching the 
quantum regime, this response can no longer be assumed 
instantaneous, therefore one must use a causal response 
function to relate the tunneling current at time t to the 
voltage at all times t' < t. It is also well known that 
the Fourier transform of a causal function is a complex 
function with real and imaginary parts related by the 
Kramers-Kronig relation. As pointed out in [17]> several 
important results are found in the quantum regime: both 
nonlinear resistive and reactive responses exist and are 
frequency dependent. The frequency dependent response is 
one of the cause of the observed mixer gain in the quantum 
regime . 
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Re[ j („)] 


= P 


oo 



— CD 


Im[j (o' )] 

u ' -u 


( 2 . 48 ) 


Im[j (u)] = I d<J (We) 


(2.49) 


Equation (2.45) completely describes the coupling 
between the quasiparticle tunneling current and the 
applied external field. As shown in (2.49) > all the 
details of the superconducting quasiparticle tunnel 
junction given in (2.46) are governed by the quasiparticle 
dc I-V characteristic, where the external driving voltage, 
V(t), enters solely through the function W(u). In the 
next section, the average current expression (2.45) will 
be used to derive the explicit quantum expressions for the 
elements of the conversion admittance matrix for the 
tunnel junction mixer. 
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2 . 4 Quantum mixer Theory 
2.4.1 Introduction 

The heterodyne receiver model described in section 
2.2 can be adopted in the quantum generalization of the 
tunnel junction mixer theory. That is, in the quantum 
theory, the nonlinear tunnel junction is still represented 
as a multifrequency multiport network which is 
characterized by the admittance matrix Y (2.9) • However, 
the elements of this matrix will now be calculated from 
the dynamic tunneling model described in section 2.3*2. 

2.4*2 Quantum Expression for the Conversion 
Admittance Matrix 

The voltage and current across a tunnel junction in 
the presence of a dc bias and local oscillator drive can 
be written as: 


v(t) = V 0 + V L 0 (t) 


(2.50) 


<I(t)> = I L0 (t) 


(2.51 ) 
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Assuming a periodic local oscillator drive, the induced 
average tunneling current (2.51) is: 


<I(t)> = I I0 (t) 


°° imu t 
= Im J e p 

m=-oo 


* 2 w L0 (nu )V/^ 0 ( (n+m) u ) 3 (nu +eV 0 /h) (2.52) 

n — —.co s' r * 


Where the phase factor (2.33): 


■/ 


-inu t 


exp[-i e/fi I dt'V L0 (t')] = ^W Tn (nto^)e P (2.53) 


n=- c 


LCr p' 


has been used. Including the small signal v . (t) and 

s 1 ^ 

i si g(t) in (2.50) and (2.51), we have: 


- v 0 + V I0 (t) 


V . (t ) 
Slg' 


(2.54) 


<I(t)> = I LQ (t) + i sig (t) 


(2.55) 
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Where the small signal v (t) and i . (t) can be 

S 1 ^ sig 

expressed as sideband frequency components: 


“ i u t 

v Bi (t) = Re I V e 11 

D EQ=— oo 


1 


-iu t i to t 


-r Z <V " + V “ > 

d m=— <» 


( 2 - 56 ) 


= Re 2 i 


00 i to t 
— m 


sig 


m=-“ 


m 


1 “ *. -iu t iu t 

= - 2 (V " + i„e m ) 


2 m=-° 


m 


(2.57) 


To the lowest order of V , the Fourier components of the 
phase factor (2.53) can be approximated by: 


W( V = z w M 

n=-» 


(no ) { 6(o p -nco p ) 
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+ 


2fiu , 
m ' 


Tv , {(u -nu -u , ) 
L m' ' p p m' 


- v , s ( u -nu +u . ) ] i 
m' v p p m ' ' J 1 


(2.58) 


the resultant current across the junction is then given 
by: 


<I(t)> = Im{W L0 (nu p )V L0 (n' u p ) 


x j(nu + eV 0 /-h)exp[-i(n-n' )u t] 


eY 

+ i, -r 2 — 2 , w i,o (no p )w io (n, “p > 

m' 2hu . nn' * * 

m 

x [j (nUp+eV 0 /ft) - j (nUp-u m +eV 0 /-h) ] 


x exp[-i(nup-n’ Up-u m , )t] 


eV* , * 

+1 — — W L0 (nu n )W 10 (n ’ u n ) 

m' 2hu , L0 P L0 P 
m 


X [d*(nUp+u m ,+eV- 0 /h) - j*(n' Up+eV 0 /-h) ] 


exp[-i ( nUp-n' u p +u m )t] 


X 
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(2.59) 

+ [higher order terms of v and v*ll 

m m 

On the other hand, the average tunneling current (2.55) 
can he written as: 



Where the small signal current-voltage relation (2.9) has 

been used. Substituting (2.52) into (2.60) and comparing 

with (2.59), the matrix elements Y . = G , + iB , can 

mm' mm' 0 mm' 

be identified as: 



X 
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- [j*(n’ u p +u m ,+eV 0 /-h) - j*(n' u p +eV 0 /-h) ] j 


and 


( 2.61 ) 


B. = Re(-i) 

nun 1 ' 


^ W T.n( nU T,) V T.n(n' “ ) 6 


2hu , nn’=- c 
m 


L0 V “' p' L0 V “ p' m-m’,n'-n 


^ { [j (nu p +eV 0 /ti) - j (nu p -u m ,+eV 0 /fi) ] 

- [j*(n' u p +u m ,+eV 0 /fi) - j*(n' u p +eV 0 /ti) ] } 


Where B mm » is the reactive mixing term which only 
appears in quantum theory and has no classical analogue. 
Replacing the matrix elements Y mm , (2.10) by (2.61), the 
performance of the tunnel junction mixer can then be 
evaluated exactly as described In the generalized 
heterodyne receiver models (2.13 - 2.25)* 
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2 . 5 SIS Mixer Noise Theory 


2.5*1 Intro duction 


The sources of noise in a practical SIS mixer are (1) 
shot noise due to local oscillator current I-^Ct) (2.51) 
tunneling through the junction barrier, (2) thermal noise 
generated from each sideband termination, and (3) Noise 
arising from the quantum fluctuations of the input 
radiation fields. Since the complete quantum mechanical 
analysis, which quantizes "both" the input radiation 
fields and the charge carriers, is not yet developed for 
tunnel junction mixer, and the noise due to the 
fluctuations of the applied radiation field is small, tiu/k 
= 5-5 K at 115 GHz, it is omitted in this analysis. 

Shot and thermal noise sources can be included in the 
mixer equivalent circuit by connecting an equivalent noise 
current source to each sideband termination as shown in 
Fig. 2-2. The equivalent IF output noise current due to 
the noise current at all the sideband frequencies can be 
determined by equation (2.19)* 


2.5-2 Shot Noise 
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In a practical SIS mixer, the mixer noise should he 
dominated by shot noise which is due to fluctuations in 
the local oscillator current I^gCt) tunneling through the 
junction barrier. Rogovin and Scalapino [50] have studied 
the current fluctuation phenomenon in tunnel junctions and 
have found that it can be simply related to the junction's 
dc I-V characteristic. Tucker [75] extented this 
calculation and constructed a complete shot noise model 
for the tunnel junction mixer. 

In Tucker's model, shot noise is represented by a 
noise generator [ I ( t ) - <I(t)>] placed in parallel with an 
ideal "noiseless" mixer. I(t) is the junction's current 
operator over a time interval *T', and <I(t)> is the 
expectation value of I(t) . As shown in Pig. 2-2, this 
noise generator can be replaced by a set of current 
generators I m (t), one at each sideband port. In that 
case: 


u+irB 
m 


I (t) « 
m 


j do’[I T (u' 


-iu't iu't 

)e + Im(-u' )e ] 


u -irB 
m 


- 2B[l T (o m )e 


“i "t iu 

■ + v-v* m ] 


(2.62) 
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where I^(u) is the Fourier component of the current 
operator I(t) defined as: 


I T (u) 



dt' iut ’ 


I(t' )e 

2 IT 


(2.63) 


The effective noise current operator at the output 
frequency Uq is given by (2.19): 


I 


eff 

T 



m ' 


om 


« ^fp( u m‘ ^ 


(2.64) 


At the output frequency uq, the mean square noise current 
within the bandwidth B is: 


T/2 

> L0 = lim 4/ dt<[l 0 (t)] 2 > (2.65) 

-T/2 


Where ^(t) can be calculated by equation (2.62). 
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Substituting equation (2.64) into (2.65), we have: 


<u n ] 2 > 


0 J LO 


B lim — <[l® ff (u Q ), lJ ff (-o n )] J .> 


T-^a> T 


0 y J +" 


= B y X- i H , 
^ , Om Om ’ mm ' 
m,m 


( 2 . 66 ) 


Where [ ], denotes an anticommutator and H , is the 

+ nun 

current correlation matrix which is defined as: 




H 


mm' 


= lim 


<[!„.<— J. I-nC ) 3^> 


(2-67) 


For a tunnel junction in the presence of a periodic local 
oscillator drive, the current correlation matrix (2.67) 
can be expressed as [73]: 


H 


mm' 


n,n' =- 


W L0^ ntJ p^ W L0^ n ' 


u ) 6 
p' m-m 


,n 


-n 


e 
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x f coth[ 8 ( eV Q + n’tiu + Hu ,)/2] 

P ® 


X 


J dc^ V 0 


+ n' Hu / e + Hu , /e) 
p m 


+ coth[e(eV 0 + nHu - Hu ,)/2] 

P ® 


x I dc (V 0 + nhu p / e " (2.68) 


where W-^nu) is the Fourier component of the phase factor 
given by equation (2.53) » and B = 1/kT. 

The equivalent USB current source I which 

s 

corresponds to the same mean square noise current <[Iq]^> 
at the output frequency can be calculated from (2.19): 



(2.69) 


and the power corresponding to such fluctuations is: 
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pshot 



4 G i I x o<\ 


<[lnl 2 > 


0 J LO 


2.5.3 Thermal Noise 


(2.70) 


The thermal noise generated from the m-th sideband 
termination can be represented by an equivalent current 
source with mean square amplitude [73]: 


<[i ] 2> 

m J 


4ft uB 

m 

[exp(fiu m /kT B )-1 ] x R 


(2.71 ) 


Where R is the resistive termination, Tp is the physical 
temperature of the termination and B is the effective 
bandwidth at output frequency. These noise components 
will be down-converted into the IF output port and are 
equivalent to a single noise current source at IF with 
mean square amplitude given by: 
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< ^ 1 Q 2 ^thermal " B m^O, 1 ' X Om 1 2 


4ft u. 


m 


[exp(ftu m /kT B )-1 ] x R 


( 2 . 72 ) 


As usual, noise generated by signal and load terminations 
are not included. 

Referred to the upper sideband, the signal power 
corresponding to these fluctuations can be calculated as 
in equation (2.70): 


pthermal _ . r . -i2. 

Ar ,, ,2 thermal 

4G 1 1 X 01 1 


(2.73) 


2.5*4 Total SIS Mixer Noise 


Prom the previous sections, the minimum detectable 
power of the SIS mixer is: 


pmin _ pShot + pthermal 
det 


(2.74) 


The sensitivity of the mixer is usually characterized 
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in terms of the equivalent input noise temperature of 
the mixer, where T.„ is defined as the temperature of the 
input termination of a noiseless but otherwise identical 
mixer, which would deliver to the IP load the same noise 
power as the actual mixer delivers when it is connected to 
a noiseless input termination. According to this 
definition, the noise temperature refers to upper sideband 
of a practical SIS mixer is given by: 


T 


M 


p min 

det 

kB 


{<[!<>] ^10 + ^ 


0 - 


thermal 


} 


4kBG.| | Aq.j | 


( 2 . 75 ) 


Here, the Rayleigh- Jeans limit has been used to define the 
equivalent temperature of the input termination from the 
minimum detectable power. 
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2.6 The 3-Frequency Approximation 

The simplest applicable heterodyne receiver model is 
the 3-port Y mixer model in which we only consider the USB 
frequency (u^), LSB frequency and IF frequency 

(uq). The local oscillator waveform across the junction 
is taken to be sinusoidal: 


V(t) = V 0 + V L0 cos(u p t) 


(2.76) 


and all higher harmonics and their sidebands are assumed 
to be shorted either by the junction capacitance or by an 
external tuning circuit. The mixer conversion admittance 
matrix (2.9) is then reduced to a complex 3 * 3 matrix. 
Assuming small IF frequency, Uq << u , and * u = u _i » 
the admittance matrix elements: 


Y ik = G ik + jB ik , i,k= -1,0,1 (2.77) 


have been explicitly evaluated by Tucker [75] as: 
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G 00 *. I_ J n (a) S„(V n +nft Ur ,/e) 


n=-» 


dV 'dc v O p' 
QV 0 


G 10 ~ G -10 


= - S nJ n(»>^- I ao (T 0 +n, 'V e ) 

0 


G 01 = G 0-1 


= -— I J n (a)[j n-1 (a) - J n+1 (a)]l dc ( V nftu p /e) 
tuj n=— 50 r 

P 


G 11 = G - 1-1 


( 2 . 78 ) 


2 flu n=-°° 
P 


I - J Ll< a >]W T 0 + “*V e) 


G 1-1 " G _H 


e <*> 


2fiUp n=-“ 


I J n (a)[J n _ 2 (a) - J n+2 («)H dc (V 0+ nhu /e) 


B io 


= “ S J n (a)[j n-1 (a) - J n+1 (a)] — I KK (V 0 +nhu p /e) 


2 n=-°° 


dV, 


B 11 - _B -i_i 
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e <*> 


2R(j n=-“ 
P 


I [Jn_l(«)-2J„(») +J n + l( a )] I KK( V 0 +n,1 “p/ e) 


B 1-1 


e <*> 


2hu n=-“ 
P 


2 [J n _2 (a)J n (a) ' 2J n-1 ( “ )J n+1 (a) 


+ J n (o)J n + 2( a) ^KK (V 0 +ntl "p /e) 


B 00 ■ B 01 ■ B 0-1 “ 0 


Where J (a) is the Bessel function of the first kind of 
n 

order n, and the argument a is given hy: 


a 



(2.79) 


I d0 (V) is the junction's dc I-V characteristic in the 
absence of local oscillator power, and 


I (V) 
KK V ' 


1 

- P 

IT 


GO 



— OD 




(2.80) 
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is the Kramers-Kronig transform of the c ( V ) and is the 

real part of the current response function (2.48). 

The dc component of the tunneling current (with the 
10 applied) is given by: 

I 0 = I J^( a )l dc (V 0 +nhu /e) (2.81) 

n=-~ * 

and the component at local oscillator frequency u is: 


Ip = - I nJ n(“) I dc (V 0 +ntltJ n /e) (2.82) 

x a n= — 00 ^ 

+ 3 f J n ( a) [ J n- 1 (a) " J n+1 (a) ^ I KK (V 0 +ntl V e) 

n=-“ 

Assume = u g and = ik are the signal and image 

* 

frequency respectively , and let Y. = Y_ = G + jB , Y . = 

1 D S S — 1 

Y i = G d - JB^ and Yq = Y^ = G^ + be the source, 
(conjugate) image and load terminations, the augmented 
impedance matrix Z' (2.16) can be expressed as: 


* As indicated by Torrey & Whitmer L 70 J and Saleh [54J, 
no loss of generality results from this assumption. 
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Z' = 


(Gn+ G s )+j ( B i i+B s ) G 10 +jB 10 G 1 _ 1 +jB 1 _ 1 


-1 


G 


01 


( G oo +G iP + ^ b i 


r 01 


G 10 “OB 10 (G 11 +G i )-j(B 11 + B.) 


(2.85) 


Once Z' is obtained, the mixer conversion loss, input 
and output impedance can be calculated from equations 
(2.22), (2.24) and (2.25) respectively. 

In the 3-port model, the explicit form of the current 

correlation matrix elements H , , which characterize the 

mm 

noise properties of the mixer, are also given by Tucker 
[73]: 

H oq = 2e f J^(a)coth[B(eV 0 +ntiu )/2]l dc (V 0 +ntiu /e) 
n=-°° ^ * 

H 10 = H -10 = H 01 = H 0-1 
2b °° 

= — I nJ^(a)coth[B(eV 0 +ntio )/2]l d (V 0 +nfio /e) 

rr M = -co * ^ 


(2.84) 
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= e f ^ J n+1 + J n-1 

n=-°° 

x coth[ B(eV 0 +n*to p )/2]l dc (V 0 +ntico p /e) 

H 1_1 = H _n 

00 

= 2e ^ J n+1 (a)J n _ 1 (a)coth[B(eV 0 +n-hu )/2] 
n=-« ^ 

x I dc (V 0 +ntltJ p /e) 

Therefore, in the 3-port model, the shot noise 
contribution to the input mixer noise temperature is given 
as : 

T ^ h0t = "ikh I* ^ 1 x oi |2 + 1 x o-i - 2 ) h i 1 

4k | X Q1 | 

+ 2(Re[x 01 ] + Ref^.,])^ + H Q0 (2.85) 

+ 2(Re[ Xq 1 ]Re[ ] + I m [ ^q-| ]l m [ 1 —1 ^ 


where 
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, m = -1 , 0, 1 


( 2 . 86 ) 


In the 3-port model, only the image termination 
generates thermal noise, which is given by: 


.thermal 


I x o-1 I ^ x ftu i G i 

4k|A Q1 | 2 exp[hu i /(kT B )-1 ] 


(2.87) 


and the total input mixer noise temperature is: 


m _ mShot ^thermal 
M “ M 


( 2 . 88 ) 


Theoretical calculations based on the formulae 
developed in this section will be used later to. compare 
with the experimental results. 


There are two examples of special interest which will 
be discussed below. 


2.6.1 Equal Signal and Image Termination 
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If the image termination is equal to the complex 

* 

conjugate of the signal source admittance, Y. = Y = 

1 s 

G - jB , the signal to IF conversion gain can be written as 
s s 

[ 17 ]: 



4 Woi 


[(G s +G 




> + (B s +B ir B i-i> ] 


( g oo +g l^ 


X d(g s ,b s ,g l ) 


( 2 . 89 ) 


where 


D ^ G s ,B s ,G lP = ^ G s “ ^- G 01 G 1 0 / ^ G 00 + G lP “ G 11^ 2 
+ * B s " [ G oi B io / ( G oo + G i^ " B ii ^ 2 
- |[& 1 _ 1 - G oi G 1 0 / ^ G 00 + G L ^ 2 

+ t B 1-1 - G 01 B 10 / < G oo + G 1^ 2 ^ ( 2 - 90) 

Note that the function D(G s ,B g ,G^) also depends on the 
junction dc I-V nonlinearity, dc bias (Vq) and pump 
voltage (V L0 ) through the admittance matrix elements Y mm , 
( 2.61 ). 
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It can be verified from equation (2.89) [17] that for 
a given tunnel junction mixer with dc bias Vq and pump 
voltage connected to a load conductance G^, it is 

possible to achieve infinite conversion gain if and only 
if D(G S ,B S ) = 0. 

We can determine the source admittance region in 
which a given tunnel junction mixer will exhibit infinite 
gain: 


If B ( G S » B S ) = equation (2.10) can be rearranged as 
follows : 


(G s - X Q ) 2 + (B s - Y 0 ) 2 . r 2 


(2.91) 


with 


r - G oi G io / ^ G oo + g l ^ 2 

+ [B 1 _ 1 - Gqi B 10 / ( G 00 + G L )]2 ^ 1/2 (2.92) 


and 


X 0 “ G 01 G 10 / ^ G 00 + G lP “ G 1 1 (2-93) 

Y 0 = G 01 B 10 / ^ G 00 + G lP “ B 1 1 


(2.94) 
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This is the equation of a circle in the source admittance 

plane with radius r and centered at (Xq, Yq). Here Xq, Yq 

and r only depend on the junction's dc I-V characteristic, 

dc bias point and pump voltage. Since in the real tuning 

structure G g > 0, therefore, infinite gain is attainable 

only on the portion of circle with G > 0. 

s 

The mixer's IF output impedance (2.25) is real and 
is : 


= [« 5 S + S ,,) 2 + ( B s + B ,,) 2 - ( G 2 _, + B 2 .,)] 

/ [G 00 * D(G s ,B s ,G l =0)] (2.95) 

with 

D(G s ,Bs,G l = 0 ) 

= C G S - (G 01 G 1 0 /G 1 0“ G 1 1 ) + ^s-( g oi b io/ g oo- b ii)^ 2 

(2.96) 

+ [ _1 — G 01 G 1 q / G 00 ^ 2 + ^■ B 1 - 1 “ G 01 B 10 // G 00 ^ ^ 
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Where is the differential resistance of the pumped dc 
I-V curve. Within a certain range of the parameters Vq, 
V^o» G g and B g , the function D(G s ,B g ,G L =0) < 0. That is, 
the output impedance of an SIS mixer can be infinite or 
even negative. This has been verified in the experiments 
(see section 6.3-1). 

For a given tunnel junction, a negative output 
impedance region can be determined in the source 
admittance plane. As shown above, the equation 
D(G ,B_ ,G t = 0) = 0 is again a circle on the source 
admittance plane with radius r and centered at (Xq, Yq) 
which is given by: 


r - ^ G i-i~ G oi G ic/ G oo-l 2 + [ B i -i “ G oi B i c/ G oo -^ 1 

(2.97) 


and 


X 0 " G 01 G 10 / G 00 G 1 1 (2.98) 

Y 0 = G 01 B 1 0 / G 00 " B 1 1 


(2.99) 
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is infinite on the circle. Within the circle, the 

output impedance is negative. Usually for a given 

junction, Xq < 0. However, for a real tuning structure, 

we are limited to the case requiring G > 0. Therefore, 

s 

for a given junction, the condition for achieving negative 
output impedance is 



2.6.2 Short-circuited Image Case 


Mixers operating in single sideband mode are of 
particular advantage for many astronomical applications 
since the noise coming from the celestial background and 
our atmosphere at the image frequency is not received, 
which makes the effective sensitivity of the system better 
for a single sideband receiver. Previous computer 
simulation [65] has shown that large conversion efficiency 
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and large dynamic range are possible for a single sideband 
SIS mixer. Techniques developed in previous sections can 
also be used to study this special case. 

Here, we assume that the image frequency wave is 
short circuited, that is, The mixer conversion 

admittance matrix (2.16) reduces to a 2 x 2 matrix and the 
conversion loss of the mixer can be expressed in the 
following simple form: 


L c = 


[ G B +Gl , -G 01 G, y (Gqq+G.) ] 2 +[B s +B 1 , -G 0) B, t j (G^+Gt) f 


4G s G L G o/ ( ' G 00 +G L^ 


( 2 . 102 ) 


This can be rearranged to give the same form as in (2.91), 
i.e. constant conversion loss circles in the source 
admittance plane with radius: 


r - {4L 2 G^gJ 1 /(G ( - )0 +G l ) 4 


+ 4lg l g 01 [g q1 g 1 0 /(g 00 +g l )-g 1 1 ]/(g 00 +g l ) { / 


(2.103) 


and 
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x 0 - 2 IG I G o/ G 00 +G L )2 + G 01 G 10 /(G 00 +G L ) " G 11 

(2.104) 


Y 0 " G 01 B 10 // ^ G 00 +G L^ “ B 1 1 

(2.105) 


Hence, in the shorted-image case, the source admittance in 
which a tunnel junction mixer with given parameters (G^, 

Vq and V LQ ) will exhibit constant conversion loss is 
constrained to lie on a circle in the source admittance 
plane with radius r and centered at (Xq, Y q). The 
condition G g > 0 limits the accessible part of the source 
admittance plane. 

In shorted-image case, the output admittance, Y ^ = 
G out + ^ B out’ ^ no lon g er rea Y and is given by: 


G 


out 


(G s^1 1 ) [W W 1 )-°oi G 1 p] + ( V B 1 1 ) [ g oq(V e i 1 )- G oi E 1 p] 

(G s +G 11 )2+(B s +B 11 )2 


(2.106) 
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B 


( G s*°1 1 ) ^ G 00 ^ B s +B 1 1 ^“°01 B 1 0 -^ B s +B 1 1 ^ ^ G 00 ^ B s +B 1 1 )^ 01 B 01 ^ 


out 


(G s +G 11 ) 2 +(B s +B 1 1 )2 


(2.107) 


For optimum operation, it is assumed that the output 
impedance is conjugately terminated by the load impedance, 
Zl = 2 QUt . High or even infinite gain is possible in this 
case. Within a certain range of the source admittance, 
the output conductance (2.106) can be zero or even 
negative. Following the same technique described above, 
equation (2.106) can be put into the form of equation of a 
circle with: 


r ^ G 00 G 01^ 2 + ( B 10 G 01 ^-i 1 / 2 
r = [ ^ J 


2G 


00 


and 


G 01 G 10 c 

x o = —7: G 11 


2G 


00 


G 01 B 1 0 

r 0 '-rr- 15 -- 5 ,, 


2G 


00 


(2.108) 


(2.109) 


(2.110) 


The additional constraint that G g > 0 implies that 



67 


r > Xq, therefore, the condition to achieve zero or even 
negative output conductance is: 


G 01 G 10 


2G 


„ , r (°oi G io ) 2 + ^oiW 2 , 1 ' 2 

- ®11 < 1 ZS 1 


00 


2G 


00 


( 2 . 111 ) 
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2.7 Theory of Arrays 

It has been shown [19] that the mixer performance of 
an N junction array should be the same as that of a single 
junction which has the same overall impedance provided 
that (a) current is in phase all along the array and (b) 
that all of the junctions of the array are identical. 

Under these assumptions, Tucker's theory for single 
junction can be applied to an N junction array by using 
the following transformation: 



SINGLE JUNCTION 

ARRAY 

Current Scale 

Kv) 

NI(NV) (2.112) 

Voltage Scale 

V 

NV 

Induced Photon Step 

■fto/e 

N hu/e 

Admittance Matrix 
Element 

Y mm ,(V,We) 

(NV,Nhu/e) 

Source Admittance 

Y 

s 

Y 

s 


Thus the conversion loss for the array will be the 
same as for the single junction with the same over-all 
impedance if the external impedance stays the same with 
the dc and LO voltage scaled by a factor N. On the other 
hand, the required LO power and the RP saturation power 
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2 

will be N larger for an N-j unction array than for a 

2 

single junction. The shot noise current < [Iq] >^q for 
an array has been shown to be identical to that of a 
single junction of the same over-all impedance [19]« 

Since the mixer conversion loss and the external impedance 
are the same in both cases, the equivalent input shot 
noise should be the same for array and single jucntion. 
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2 . 8 Computer Modeling 

As shown in previous sections, in the quantum mixer 
theory, the mixer performance is completely determined by 
the junction's dc I-V characteristic and the embedding 
admittances at all relevant frequencies. Based on this, a 
computer program was developed to determine the mixer 
performance. 

In these calculations, the junction capacitance is 
treated as an external circuit element and is combined 
with the embedding admittace. The 3-frequency 
approximation was used in the simulation. This is 
equivalent to assuming that all the higher harmonics and 
sidebands above are short-circuited by the 

junction's capacitance. This is a reasonable 
representation of our experiments since the junctions' 
uR^C product is usually larger than 7* 

In order to compare the theoretical predictions with 
experimental results, the experimental dc I-V curve of the 
junction was used in the simulation. Furthermore, the 
embedding admittance, seen by the SIS junction at signal, 
image and L0 frequencies, is a function of the backshort 
position and is determined by direct admittance 
measurements on a 40X scale model of the mixer block. 
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Given the input pump power from the experiment and 
the embedding admittance at the 10 frequency from scale 
model measurements, the pump voltage impressed across 
the junction can be determined from the following 
equation: 


P L0 = ~T~ [(V L0 G 10 +I U > 2+ < V I0 B I0 +I - u f 2 - 11 ^ 

8G L0 p P 

Where I (I ) is the real (imaginary) part of the 
°P ” U P 

tunneling current at the local oscillator frequency 
(2.82). Since the argument of the Bessel function is V^q 
dependent, equation (2.11 3 ) can be solved iteratively. 

With all these experimental data input, this program 
can predict the mixer conversion loss, equivalent input 
noise temperature, input and output impedance of the 
mixer. Experimental results are compared with those 
calculations in chapter 6. 

2.8.1 Outline of the Computing Method 

The calculation procedures can be roughly divided 

into three parts: (1) given the junction's dc I-V curve, 

calculate the Kramers-Kronig transform, (2) given I dc » 

^KK’ ^ as an< ^ P um P voltage, calculate the mixer 

conversion admittance matrix elements Y . , with m,m’ = 

mm 
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-1 , 0, 1, (3) given the source, image and load admittances 
together with Y calculate the mixer conversion loss, 
mixer noise temperature, input and output impedance. 

The experimental dc I-V curve was fed in at constant 
voltage increments with the incremental voltage equal to 
one tenth of the photon step, i.e. V = Hu^/ ( 1 0e) . For 
voltage between 2A/e - Hu /(2e) and 2A/e + Ho /(2e), 

& Jr 

exponential interpolation was used between the input 
points, whereas a third order polynomial was used outside 
this voltage range. Given dc I-V curve, the next task is 
to evaluate it's Kramers-Kronig transfrom: 

1 f dV' 

V v > = ; p c< v ’> ( 2 - 80 > 


This integral is divergent. However, 


^KK can 


shifted by an arbitrary linear term (corresponding to 


adding a fixed conductance across the junction) without 


affecting any observed quantities such as the tunneling 
current [80] and it is also known that above the gap 


voltage, the dc I-V curve should exponentially approach a 
straight line with a slope equal to junction's normal 


conductance G^. We can therefore subtract the ohmic 
portion from this integral and modify equation (2.80) as 
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follows [75]: 






(V*) - G n V'] 


(2.114) 


Since 1^ c ( V ) exponentially approaches the normal 
conductance line, I = G^V, let V^ AX be the voltage where: 


I dc^ V MAX^ “ G N V MAX " G » e *° 


(2.115) 


then equation (2.114) can be written as: 


i av 


i KK (v) = - 


I /• U V 

- p / — 

IT / V'-1 

_V “ 

MAX 


CW 7 ' 5 - g n v '] 


+ correction term 


(2.116) 


Usually, by choosing V MAX = 3Vg a p, the correction term in 
(2.116) is small enough to be discarded. Notice that at 
V 1 = V, the integrand in equation (2.116) diverges. To 
remove this singularity, (2.116) can be rewritten as 
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follows : 


i M i x av ’ 

W T > ■ ; I — l [I ao< v ')- G K v '] - 

—xrJ 


MAX 

i r dV 1 

Tf / V » 
-V J 


[Ih,(V) - G_V] 


- + 4l ) t v > 


(2.117) 


The second integal can be evaluated as follows: 


4k = -[i dc (v) - 


- VI 


V MAX + V ' 


and the integrand in the first integral is no longer 
singular at V* = V: 


v-v V'-V 


•IClrt^CV') - G„V'] - [I„„(V) - G„V]| 
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= lim 
V' +V 


W v ' 


) - 


1 dc< V > 


V'-V 


- lim 
V' -V 


V 


- V 


V'-V 


dI dc (V,) 


dV' 


- G 


V'=V 


N 


(2.119) 


which is the difference between the junction's real 
differential conductance and the junction's normal 
conductance and can be calculated from the junction's dc 
I-V curve. The first integral in (2.119) can therefore be 
evaluated by numerical integration. The resultant 
Kramers-Kronig transform 1^ used in the simulation is: 


1 

J KK (V) = - G M V]*n 


dc w ' U N' 


V MAX " V 


V MAX + V 


( 2 . 120 ) 


i iV 


+ - f — 

TT / V'-V 

_v* / 

V MAX 


{[ I dc (V ’ ) - V] - ^^(V) - G W V]J 




dc w/ 


Assuming an ideal SIS I-V characteristic [22,80] and 
Vjvjax = 5^gap» calculated Kramers-Kronig transform I KK 
using (2.120) is found to be within 1# of the analytical 
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results over most of the region. About 20^ error was 
found at the gap voltage point. However, since in 
experiments, the bias points of the best performance are 
always in the middle of one of the photon steps, only the 
points with V = 2A/e - tiu /(2e) ± nHu /e, n = 0, ±1, ±2, 

.r' ^ 

...., were used in calculations. Therefore this error 
won't cause any problem in the whole analysis. 


Given I„„(V), it is straight forward to calculate the 


admittance matrix elements G.. and B. . and current 
correlation matrix elements H. . from (2.78) and (2.84). 
Note that since higher order Bessel function approach zero 


for small argument, only the most important 15 terms, i.e. 
n = -7, ...., 7, in the infinite series (2.78) and (2.84) are 


counted in this calcualtion. 


Having Y.. and H. . for a particular pump amplitude 

1 J 1 J 

and bias voltage, the mixer properties can be evaluated 
for any source and load admittance. 


77 


CHAPTER J. JUNCTION DESIGN AND FABRICATION 


3 • 1 Introduction 


For practical mixer application, superconducting 
tunnel junctions should have the following properties: 
long term stability, good thermal recyclability, high 
critical temperature (which will simplify cryogenic 
cooling), high mechanical strength, good adhesion between 
metal and substrate, and low leakage current below the gap 
voltage, I-V curve as close as possible to the ideal case. 
Lead alloy junctions and niobium junctions, which have 
been shown to fullfill the above requirements, are used in 
these SIS mixer experiments. The design and fabrication 
of the junctions used in this work will be described in 
this chapter. 


3.2 Optimum Junction Parameters for Mixer Applications 


There are several parameters of the superconducting 
tunnel junction which have to be optimized to achieve the 
best mixer performance. The two most fundamental 
electrical properties of the junction are the normal state 
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resistance R^ and the inherent junction capacitance Cj. 
However, with given junction materials and fabrication 
process, the only two parameters which can be directly 
controlled in fabricating the junction are the junction 
area A, and the maximum dc Josephson current density j^. 
The focus of this section is to determine the optimum 
fabrication parameters of SIS junctions for mixer 
applications . 


3.2.1 Optimum uR^Cj Product 


In an SIS mixer, the nonlinear effect of the 
quasiparticle tunneling current is used in the mixing 
process. The inherent junction capacitance, being in 
parallel with the intrinsic junction, tends to 
short-circuit the desired junction nonlinearity and thus 
to prevent the full signal power from being coupled into 
it. It is possible in principle to tune out the junction 
capacitance at the signal frequency, thereby avoiding this 
power loss. The junction's uR^Cj product , which is the 


* Actually, the mixer input resistance R. , not R„, 

should be used here. But since they are iBmparable in 

most cases (at least of the same order), R„ is used 

instead of R. in most of the literature, 
in 
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ratio of the junction's normal state resistance to its 
capacitive reactance, is commonly used to measure the 
effect of the junction capacitance on the RF power 
coupling. 

Although the optimum value of uR^Cj is still an open 
question, previous published experiments have shown that 
small values, i.e. = 1, tend to give poor mixer 

conversion [14»15]» whereas good mixer performances were 
obtained with uR^Cj <r 10 [33,61]. It has also been 
argued that uR^Cj > 4 is required for mixer's conversion 
to reach the values predicted by Tucker's theory assuming 
the three-frequency approximation [53]* Thus we chose 
uR^Cj = 6 in designing the Princeton - NBS junction used 
in this experiment. This value should permit operation 
close to the three-frequency condition, while not being 
too difficult to tune out. 


3.2.2 Normal State Resistance 


In order to achieve optimum mixer performance, the 
junction's resistance should be chosen for reasonable 
input and output impedance matchings. Previous 
experiments and theoretical analysis have shown that the 
input resistance R^ n and the output resistance R ou ^. 


of an 
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optimized mixer are comparable to its normal state 
resistance. Since the output of the mixer is usually 
coupled to a 50 ohms IF amplifier and the guide impedance 
of the quarter hight WR-10 waveguide used in our 
experiment is about 111 ohms, therefore, the desired 
junction normal state resistance is chosen to fall in the 
range from 50 to 100 ohms. 


5.2.5 Junction Capacitance 


With the given uR^Cj product and junction normal 
state resistance, the desired junction capacitance at the 
operating frequency can then be determined. This 
capacitance must be resonanted out at the signal frequency 
for efficient RF coupling. In principle, this can always 
be achieved by using an external tuning circuit since this 
is a constant capacitance and the RF resistance in the 
superconducting electrodes is very small. By choosing 
uRjjCj = 6, the desired junction capacitance at 115 GHz is: 


C J = 


2 it X 115 X 1 0^R 


N 


81 


= 83 fF with R^ = 100 ohms 

(3.1) 

= 1 66 f F with R^ = 50 ohms 

3.2.4 Junction Area 

For a tunnel junction with barrier thickness d and 
relative dielectric constant e, the specific capacitance 
is simply given by the parallel plate capacitance formula 



(3.2) 


However, since the thickness, composition and 
structure of the oxide barrier are difficult to 
characterize, experimental data are used to determine the 
junction area. According to Magerlein [37], the specific 
capacitance for the Pb-In-Au/oxide/Pb-Bi tunnel junction 
fabricated by standard IBM process is: 


1/Cg = 0.35 ± 0.040 - (0.016 ± 0.006)- * *n(j) (3-3) 
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for junctions on the films containing 13*6 at# indium. 

2 2 

Here Cg is in pF/cm and j is in A/cm . For critical 

P 

current densities between 500 and 1500 A/cm , Cg is 

p 

approximately equal to 4-2 pF/cm . The desired junction 
area is then given by: 


A 



(3.4) 


3*2.5 Critical Current Density 


Theoretically, for a given normal state tunnel 
resistance R-, the maximum dc Josephson current can be 
obtained from Ambegaokar - Earatoff relation [ 1 ] : 


I 


C 


TT a (T) 1 a (T ) 

tanh 

2 e R n 2k T 


(3-5) 


Where a(T) is the superconducting energy gap at 
temperature T. This yields Ig = 2.15 mV/R^. for lead alloy 
junction with A(T)/e = 1.42 mV and T = 4.2 K. In 
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practice, the experimental value is 1^ = 1.6 mV/R^ for NBS 
process. Therefore, the desired junction current density 
is given by: 

1.6 x 1 0“ 3 

j = (3.6) 

R N x A 

P ? 

Here area A is in cm , R^ is in ohm and is in A/cm . 



3.2.6 Series Array 


As discussed in section 2.7, the performance of an 
N-j unction array is the same as that of a single junction 
provided that both have the same product. From 

equations (3.4) and (3*5), it can be shown that uRjjCj is 
independent of the area of the junction and depends only 
on the maximum dc Josephson current density. Therefore, 
except for the junction area, the parameters obtained in 
previous sections can be adopted in designing an 
N-junction array. The area of the individual junction 
must be scaled by a factor of N. Table 3-1 summarizes the 
design parameters for the arrays used in this experiment. 


Table 3-1 • Summary of the design parameters of the 
Princeton - NBS lead alloy junctions. 


PARAMETERS 

UNIT 

A 

B 

C 

D 

NO. OF JUNCTION 

— 

2 

2 

4 

4 

uRjjCj PRODUCT (115GHz) 

— 

6 

6 

6 

6 

CURRENT DENSITY 

■ 

A/ cm 2 

800 

800 

800 

800 

ARRAY RESISTANCE 

ohm 

50 

100 

50 

100 

ARRAY CAPACITANCE 

fF 

166 

83 

166 

83 

AREA OF 

INDIVIDUAL JUNCTION 

im ^ 

8 

4 

16 

8 














3*3 Junction Fabrication 


Junctions fabricated by different laboratories using 
similar processes were tested in this experiment. The 
14-element series-connected arrays of lead alloy tunnel 
junctions supplied by IBM were fabricated on the 0.010" x 
0.005" * 0.010" thick silicon chips by the standard IBM 
photolithographic process [28]. In this section, we describe 
the similar process used at the National Bureau of Standards 
in Boulder where junctions were fabricated to our design by P. 
Timbie of the Physics Department, Princeton University. 

3.3*1 Princeton - NBS Junctions 


The Pb-In(12 wt$)-Au(4 wt#)/oxide/Pb-Bi (29 wt$) tunnel 
junctions were made on 0.010" thick 2" diameter fused quartz 
wafer (Mark Optics, type GE 124)* The superconducting 
Pb-alloy layers and the insulating SiO films were 
vacuum-deposited on the quartz wafer and patterned using 
optical lithography and the photoresist stencial lift-off 
method. Tunnel barriers are formed on the base electrode 
films by thermal and sputter oxidation. A wafer contains 
thirty-six 0.246" x 0.246" chips. Each chip is identical in 
design and contains 4 different types of. arrays (N = 1 , 2, 4> 
8) on it. Pig. 3-1 shows a 0.010" x 0.005" 2-junction array. 
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The base electrode is formed by sequentially depositing 

films of Pb (1000 A), Au (45 A), Pb (600 A) and In (350 A) on 

_ 0 

the quartz wafer. It is then covered with a 3000 A thick SiO 
insulation layer. Prior to the SiO deposition, a photoresist 
stencil was prepared over the base electrode. The junction 
windows are opened in the SiO by the lift-off process to 
define the junction. The oxide barrier was formed on the 
base electrode by sputter-etching in an RF-oxygen plasma 

0 

(P s 16.7 u ) . After the tunnel barrier was formed, a 4250 A 
thick Pb-Bi(29 wt$) counter electrode and a 1000 A SiO 
protection layer were deposited. An additional 1 um thick SiO 
protection layer was put down to cover the whole of each 
device except for the bonding pads. Finally, the entire wafer 
was covered by 1 urn thick photoresist, for protection, and the 
fabrication process is completed. 

All the single-junctions we obtained were short 
circuited, probably due to the indium contamination of the 
Pb-Bi electrode during the fabrication process. None of the 
8-junction arrays has been tested due to the wrong spacing 
between the two bonding pads. Therefore only 2- and 
4-junction arrays have been tested. In general, most of the 
2- and 4-junction devices show very similar I-V 
characteristic. However, arrays with an anomalous junction, 
which has a totally different critical current density from 
the rest of the junctions in the same array, are found in some 
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chips. The RF performance of these anomalous arrays is in 
general inferior to that of the normal ones, there is no 
attempt in this work to look into this problem. 

3 • 4 Junction Storage 


It has been reported that the best storage condition for 
the lead alloy junctions is to keep them in liquid nitrogen 
[5]« However, this is impractical for our application. In 
our laboratory, the junctions are stored in a dry nitrogen 
atmosphere at a temperature of 0 C. Chips covered by 
photoresist have remained essentially unchanged for more than 
two years. Two arrays with photoresist removed did fail, and 
exhibited an S-I-N type I-V characteristic after 6 months from 
the initial cool down to liquid helium temperature. But in 
general, the junctions' thermal recyclability is good. In 
particular, one 2-junction array has been cycled more than 15 
times between room temperature and liquid helium temperature 
over a period of 14 months with no discernible change in 


behavior . 
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CHAPTER 4 . MIXER BLOCK DESIGN AND 
CHARACTERIZATION 


4.1 Introduction 


As we have seen in chapter 2, the performance of an 
SIS mixer is completely determined by the junction's dc 
I-V characteristic, input LO power, bias voltage, and 
embedding admittance seen by the junction at the local 
oscillator and all relevant harmonics and sideband 
frequencies. This embedding admittance is mainly governed 
by the waveguide and microstrip circuits immediately 
surrounding the junction, i.e. the "mixer block". Thus, 
one of the main tasks in developing a high frequency SIS 
mixer with good performance is to design an external 
microwave tuning circuit (mixer block) which can provide 
the junction with the optimum embedding admittance. 

In this chapter, we describe the mixers used in this 
work. Techniques used to characterize the mixer mount 
will be described, and a microwave equivalent circuit will 
be derived which includes the effects of conductor loss. 
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4*2 Mixer Block Design and Construction 


Two different blocks were used in this work. The 
first design. Type J-J, was used in peliminary tests to 
verify Tucker's theory. The second one, Type-D, which has 
greater tunability was used in subsequent experiments to 
achieve better mixer performance. 


4.2.1 Type J-J Mixer 


This mixer mount, described by Y. Taur and A.R. Kerr 
[66], was originally designed for point-contact Josephson 
junction mixers. The mount, shown in Pig. 4-1 (a), 
consists of two major parts, the waveguide transformer and 
the main body. The transformer reduces the full height 
WR-10 waveguide (0.100" x 0.050") to quarter height 
(0.100" x 0.0125") in a series of quarter-wave steps. It 
is electroformed copper. The main body of the mixer is a 
brass block split along the wide wall of the waveguide 
with a 0.019" wide x 0.017" deep slot which contains a 
suspended quarts substrate. Metal choke patterns are 
fabricated on the quartz substrate by photolithographic 
techniques. The metallization is 2.5 um thick gold 
over a thin (80 A) chromium layer for better adhesion. 
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As shown in Fig. 4-1 (b), there are two separate metalized 
patterns on the substrate, the upper (lower) one consists 
of 5 (3) quarter wavelength stripline sections of, 
alternately, high and low characteristic impedance, and an 
arrowhead-shaped section. These RF choke patterns are 
designed using large scale- models as described later in 
this chapter: the arrowhead-shaped elements are trimmed on 
the scale model until a proper reactive impedance for 
optimum coupling to the junction is measured. This final 
shape is used in the actual mixer. 

Two 0.100" long, 0.001 "-diameter gold wires are 
soldered on the upper end of the substrate before the 
assembly for later connection to the dc/lF connector. The 
SIS chip, described in section 3-3, is mounted on the 
substrate by cold welding with indium. The assembled 
chip-substrate unit is then transferred into the mixer 
mount. As shown in Fig. 4-1 (c), it lies sideways in the 
substrate slot of the block and is centered in the slot by 
a spring-loaded clamping wire. The same spring wire also 
presses the lower part of the substrate against an indium 
coated shoulder in the block, to provide dc and IF return 
paths, and an RF ground. An adjustable contacting type 
backshort behind the chip is used as the sole tuning 


element . 


4.2.2 Type-D Mount 


Figure 4-2 is schematic drawing of the type-D mixer. 
The gold plated brass mixer block consists of two parallel 
quarter height WR-10 waveguides. It is fabricated in two 
halves with the main and side-arm waveguides split along 
the center of their broad walls. Since at 110 - 118 GHz 
we are operating somewhat beyond the high end of the 
normal operating band of WR-10 waveguide (75 - 110 GHz), 
the actual width of the waveguide was reduced from 100 
mils to 96 mils to prevent unwanted excitation of higher 
waveguide modes. As shown in Fig. 4-2(a), transistions 
between the main waveguide and the input full height 
waveguide are accomplished by a broad-band channel 
waveguide transformer [59]- A 0.020" x 0.010" substrate 
channel supports a suspended substrate stripline circuit 
between the two waveguides. The substrate, shown in Fig. 
4-2(b), is 0.023" wide x 0.003" thick fused quartz, 
metalized with stripline circuit patterns which consist of 
an RF choke structure and a 0.059" long stripline of 50 
ohm characteristic impedance. The RF choke is a 7 element 
quarter-wave design which is used to pass the IF and dc 
bias to the junction while preventing RF power leakage. 

The side-arm waveguide is coupled to the SIS chip in the 
main waveguide through the 50-ohm stripline, thereby 
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providing a tunable reactance in series with the chip. 

Two contacting finger type moveable short circuits, shown 
in Pig. 4-2(a), are used as tuning elements. The side-arm 
waveguide is blocked by a fixed indium short circuit a 
quarter wavelength away from the substrate channel; this 
reflects an open circuit at the stripline-waveguide 
transistion. The lower half block has 5 mil surface 
relief over the shaded area, shown in figure 4-2(a), to 
ensure a good contact along the joint line and reduce loss 
due to leakage. Gold wires, 0.005" in diameter, are 
placed in the gasket groves around the SIS chip area to 
prevent microwave power leakage through the interface. 
Compared with the conventional mixer design with only one 
tuning element, this design has much greater tunability at 
a given frequency, but the tuning is very frequency 
dependent . 


4.2.5 Mixer Assembly 


The mixer assembly is quite simple. Referring to 
Pig. 4-2(a), the two 0.100" long 1 mil diameter gold wires 
are first soldered to the substrate on the choke pattern 
end, a 0.050" x 0.009" * 0.001" gold grounding tab is then 
bonded on the other end of the substrate. Two tiny pieces 
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of indium are melted onto the middle of the substrate 
where the SIS chip is to be bonded and then scraped with a 
blade to expose a fresh, level surface. The SIS chip, 10 
x 5 x 10 mil thick, is then cold welded to the indium pads 
on the substrate using a pressure mounting jig. The 
complete chip-substrate assembly is transferred into the 
mixer mount. The substrate is located such that the SIS 
chip sits in the main waveguide while the lower end of the 
substrate extends across the side-arm waveguide with the 
grounding tab resting on the split surface of the block. 
The two gold connecting wires are then attached into the 
indium coated center pin of the SMA connector for junction 
dc bias and IF output. The upper half block is then 
bolted in place. When the block is assembled, the 
grounding pad on the substrate is squashed into a 0.020" x 
0.005" x 0.004"-deep slot in the upper half of the block. 
This slot is filled with indium and covered with a 1 mil 
gold foil to ensure a rigid support for the substrate and 
a good electric ground contact. Finally, the backshorts 
are slipped into the waveguides and the mixer is ready for 


use. 
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4 . 3 Mixer Mount Characterization - Scale Model 
Measurements 


4.3*1 Introduction 


Due to the complexity of the mixer geometry, 
theoretical characterization of the mixer would be 
extremely difficult. Direct impedance measurements are 
very hard to perform accurately above about 10 GHz. Since 
the physical dimensions of the mixer mount can be scaled 
inversely with frequency without appreciably affecting its 
electric properties (except for the conduction loss), it 
is much easier to characterize the mixer mount by using a 
large-scale model at a frequency much lower than the 
normal operating frequency. The loss in the mount, which 
is due to the skin effect in the conductors, can either be 
modeled by choosing the appropriate construction material 
for the scaled mount or by theoretical estimation. Scale 
modeling is a valuable technique commonly used in 
designing millimeter-wave components. See references [25] 
and [58]. 


4.3*2 The 40X Scale Model 
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Two aluminum mixer models were made, one for the type 
J-J mount, and the other for the type-D mount. In this 
section, we shall thoroughly discuss the characterization 
of the type-D mount. Since the same technique is used to 
characterize the type J-J block, we shall just present, i 
Pig. 4-17, the measurement results on type J-J model 
without going into the details. 

For the mixer model, a scaling factor of 40 was 

chosen so that the operating frequency scales from 116 — 

112 GHz to 2.9 - 2.8 GHz. Since the junction's 

capacitance is large enough ( uR^C ^ 7 to 10) effectively 

to short circuit all the higher harmonics, only the 

embedding admittance in the fundamental frequency band 

needs to be measured. A cross-section of the scale model 

* 

is shown in Fig. 4-3* The substrate is made of stycast 

which has the same dielectric constant (e r = 5*8) as fused 

** 

quartz. Copper adhesive tape on the substrate surface 
was cut to form the required circuit pattern. As shown in 
figure 4-3> a 50 ohm 0.085" diameter coaxial cable is 
buried within the stycast substrate in contact with the 
copper tape, with one end connected to the network 


Emerson & Cuming Inc. Canton. MA., Stycast Hik. 

** 3-M Co. St. Paul, MINN. Electrical tape No. 1245* 
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analyzer. The other end is at the position of the SIS 
chip in the real mixer, and hence enables us to measure 
the embedding admittance normally seen by the junction. 
Some small corrections must be made, as described later. 

4*5*5 Measurement Setup 


Figure 4-4 is the block diagram of the measurement 
setup. The test signal generated by the sweep oscillator 
(HP-8690B) is coupled to the power spliter (HP-1 61 OA), two 
identical 10 dB attenuators (HP-8492A) and two 
high-directivity directional couplers (Narda-5295) • One 
directional coupler is connected to the scale model, the 
other connected to a short circuited reference coaxial 
cable identical in length to the cable into the model. A 
Hewlett Packard network analyzer (HP841 1 A/841 OA ) compares 
the reflected signal from the reference cable with the 
reflected wave from the mixer model and displays the 
embedding admittance seen by the end of the coaxial probe 
(corresponding to the SIS chip in the real mixer) directly 
on the polar display (HP-841 4A). The results can also be 
recorded on a Smith chart by an X-Y recorder (HP-7034A). 

4*5*4 Corrections to the Measured Embedding 


Admittance 
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Since the coaxial probe does not resemble the SIS 
junction (which is a planar device) in the chip region, 
several small corrections have to be made to transform the 
measured embedding admittance (on the scale model) into 
the embedding admittance seen by the junction in the 
actual mixer. This can be accomplished by characterizing 
the SIS chip and the test probe as follows. 

Figure 4-5 is an enlarged view of the SIS chip in the 
actual mixer and the corresponding region in the mixer 
model. Since the physical dimensions of the SIS chip in 
the actual mixer and the test probe in the scale model are 
much smaller than the guide wavelength, 121.6 mil at 115 
GHz and 4.86" on the 40X scale model, they can be 
characterized by a number of lumped elements between the 
reference terminal planes A-A', B-B' shown in Figs. 4-5(a) 
and (b). The test probe can be represented by the 
inductance of the center conductor L' and the capacitance 
C' at the end of the coaxial cable, while the SIS chip can 
be characterized by the capacitance C and the equivalent 

s' 

inductances Lj and Ljj which represent the magnetic energy 
stored in regions I and II shown in Fig. 4-5(a). Three 
large scale models were constructed to measure these 
quantities, they are: 1) a 300X model of the actual SIS 
chip (similar to Fig. 4-5(a)), 2) a 7.5X model of the test 
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probe in the 40X scale model (40 x 7*5 = 300) (similar to 
Pig. 4-5(b)), and 3) the reference model which is similar 
to the 300X model of the SIS chip except that, as shown in 
Pig. 4-5(c), the conductors between the terminal surface 
AA'-BB' are removed. 

The capacitance measurements were performed with a 
Boonton Capacitance bridge at 1 MHz. is the difference 

in the measured fringing capacitance betwen the 300X model 
of the SIS chip and the reference model, while is the 
difference between the 7.5X model of test probe and the 
reference model. 

Figure 4-6 is a block diagram of the setup used to 
measure the equivalent inductances , Bjj. A low 
frequency (10 - 20 MHz) signal generated by the oscillator 
is transmitted through the 50 ohm system to the power 
meter. The unknown inductance (i.e. model of the chip) is 
then connected in shunt with the 50-ohm line. By 
measuring the change in the transmitted power the unknown 
inductance can be determined fairly accurately. A 
detailed description of this method is given in Appendix 
I. 


The measured values of these lumped equivalent 
circuit elements are summarized in table 4-1 • The 
corrected embedding admittances can be computed from the 



99 


circuit illustrated in Fig. 4-7* The procedure can be 
regarded as de-embedding the end of the coaxial measuring 
probe (-L' , -C ' ) , then re-embedding it (L^, V hi) in 
the actual SIS chip geometry. 


Table 4-1 Measured Values of the Lumped Elements 



Test 

Probe 

SIS Chip 


L' (nH) 

C f (pF) 

bj (nH) 

(nH) 

C p (pF) 

300X 

Model 

7*81 

0.462 

5*44 

0.94 

0.576 

■5SII 

0.026 

0.00154 

0.018 

0.0031 

0.00192 


4.3*5 Available Embedding Admittance: Measured 
Results 


Embedding admittances (junction capacitance excluded) 
measured on the mixer model at 2.875 GHz (equivalent to 
115 GHz in the actual mixer) are plotted as a function of 
the two backshort positions in Fig. 4-8. The results are 
shown on a Smith chart (with both conductance and 
susceptance normalized to 1/(50 ohm)). It is found that 
except for a small forbidden region (shown on the Smith 
chart), all the other embedding admittances can be 
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achieved by varying the two backshort settings. 

To obtain the embedding admittance seen by the 
intrinsic junction, the junctions capacitive susceptance 
must be added in parallel with the measured embedding 
admittance, that is: 

Y emb = Y emb + 3 “ C < 4 -’> 

This has the effect of moving each point in the forbidden 
region clockwise along a constant conductance line. The 
resultant forbidden region as seen by the intrinsic 
junction (which excludes the junction capacitance) is 
illustrated in Pig. 4-9* Therefore, almost any desired 
embedding admittance can be presented to the intrinsic 
junction in this mount. 
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4.4 Theoretical Characterization of the Mixer Mount 


4*4.1 Introduction 


Although the junction's embedding admittance can be 
determined from the scale model measurements described in 
the previous section. This technique does not simulate 
the loss in the mixer mount and is very time-consuming. 

It would be very helpful to construct the equivalent 
circuit of the mixer mount, including loss, and employ it 
to predict the junction's embedding admittance as a 
function of the two backshort positions. 

4*4.2 Microwave Equivalent Circuit 


An approximate microwave equivalent circuit of the 
type-D mount for the dominant TE^ q rectangular waveguide 
mode is illustrated in Fig. 4-10. The transistion between 
the junction and the single-mode waveguide is represented 
by a simple T-network. The elements in this network, 
which represent the electric and magnetic energy stored in 
the various parts of the mixer mount, are in general 
complex. By assuming a loss-less network, these elements 
can be identified as the gap capacitance G and the lead 
inductances & Lg as shown in Fig. 4-11. The junction's 





102 


embedding admittance can be calculated if the precise 
values of these lumped elements are known. The same 40X 
scale model described in the previous section, with a 
small modification of the printed circuit on the 
substrate, can be used to measure these elements. Details 
of the measurements are described in Appendix II. The 
measured values of these parameters are summarized in 
Table 4-2. 


Table 4-2 Measured Values of Equivalent Circuit 
Parameters 


L 1 (nH) 

L 2 (nH) 

C(fP) 

0.111±0. 019 

0. 039±0. 01 7 

6. 9±1 .01 


A simplified equivalent circuit with measured lumped 
elements is shown in Pig. 4-12. The main backshort is 
represented by the variable shunt susceptance B-g/g while 
the side arm backshort acts as a variable reactance 
in series with the junction in the circuit. The backshort 
impedance is given by : 

jXg/g = j Zg x tan(2irS/ A ) , (4.2) 

where S is the distance between the backshort and the 
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junction, is the guide-wavelength, Zq is the waveguide 
impedance of the q mode given by: 

Z G = (h/a)2407r/[l-(f/f c ) 2 ] 1 / 2 (4-3) 

Where a and b are the waveguide dimensions and f is its 
cut-off frequency. The SIS junction is represented by its 
intrinsic junction shunted by junction's capacitance. 

In Pig. 4-1 3> we superimpose the forbidden embedding 
admittance region, derived from the equivalent circuit, on 
the experimental results measured on this model. The 
slight difference between the measured and calculated 
forbidden region is due to the errors in the measured 
values of 1^, and C, resulting from the errors in the 
mechanical tolerances in the model and the errors within 
the network analyzer. It is clear from this comparison 
that the equivalent circuit is adequate to represent the 
real mixer mount. 

4.4*3 Loss in the Mixer Mount 


The maj or loss components in the mixer mount are the 
sliding backshorts, the waveguide, and the 50 ohm 
suspended stripline which couples main and side-arm 
waveguides. A movable contacting backshort usually 
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presents a small resistance between the contacting fingers 
and the waveguide walls. For instance, a backshort 
resistance as high as 7*6 ohm (an extreme case) has been 
reported in a 115 GHz Schottky diode mixer using 
backshorts similar to those we used in the type-D mount 
[23]* In order to reduce this loss, special care has been 
taken during the mixer assembly: mixer mount and 
backshorts are cleaned thoroughly; the backshort' s fingers 
are adjusted to ensure a good contact with the waveguide 
walls. The backshort resistance is estimated to be 0.1 
ohm in the type-D mount. 

The waveguide attenuation due to the normal skin 
depth effect can be characterized theoretically [10]. At 
room temperature, the theoretical loss of the gold plated 
quarter-height WR-10 waveguide is 0.245 dB/inch at 114 
GHz. At liquid helium temperature, this attenuation is 
estimated to be 0.076dB/inch . However, due to roughness 


* Since the anomalous skin effect occurs roughly near 
1 30K for gold at 114 GHz [38], the extreme anomalous 
surface resistance R m = ( ir^f^ii^/ek) ^ ^ is employed in 
this estimate. Here 6 = 7*2 is a constant, f is the 
operating frequency, n is the permeability and k is the 
ratio of the dc conductivity to mean free path a/x which 
is approximately equal to 8.4 * 10^ (ohm m^) - ^ for gold. 
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of the surface of the waveguide, the actual loss may he 
higher than this values. 

The major loss components in the 50 ohm stripline 
include the conductor loss and dielectric loss [30]. The 
attenuation due to the dielectric loss is about 0.05 
dB/inch at 100 GHz for the fused quartz substrate with 
dielectric constant e r = 3*8 and loss tangent tan 6 = 
0.0001 [30]. Regarding the conductor loss, computation 
using the analysis described by Schneider [55] yields a 
conductor loss of 1 . 96 dB/inch for room temperature 
stripline at 1 1 4 GHz and 0.65 dB/inch when cooled down to 
4 K. Since these are very rough estimates and other 
losses such as the loss due to roughness of the surface of 
the substrate channel are not included, the actual loss is 
certainly much higher than the calculated values. 

Therefore an estimated total loss of 1.2 dB/inch is 
adopted for the stripline at 4 K. 

4.4.4 Embedding Admittance Calculated from 
Equivalent Circuit 


Given the equivalent circuit of the mixer mount and 
the loss in the circuit, the embedding admittance seen by 
the junction at signal, image and pump frequencies can be 
calculated as a function of the two backshort settings. 
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This was accomplished by using a network analysis program 
[34] which computes the following quantities: a) input 
impedance (admittance) , b) input reflection coefficient 
| p | and Arg (p), c) input return loss, d) transmission 
loss. The results are printed out and displayed on the 
Smith chart directly. 

Some of the possible tuning conditions for the 
short-circuited image case, as simulated by the program, 
are presented in Fig. 4-14 to 4-16 : the embedding 
admittance seen by the junction at signal, image and 10 
frequencies are plotted on the Smith chart with main 
backshort position as parameter. Also shown is the 
transmission loss of the circuit at each possible 
embedding admittance. A transmission loss of 1-3 dB is 
usually found in the vicinity of the best operating point. 
As discussed in chapter 6, the optimum theoretical mixer 
performance is degraded by this loss. 

From the simulation, it is found that tuning is very 
sensitive to the backshort settings; an accuracy better 
than 0.5 mil in backshort position is required to achieve 
optimum results. 
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CHAPTER 5. EXPERIMENTAL SET-UP FOR SIS MIXER MEASUREMENT 


5 • 1 Introduction 


In order to evaluate the performance of the SIS 
mixer, a carefully calibrated microwave cryogenic mixer 
testing set-up was used in the experiment. Fig. 5-1 is a 
block diagram of the experimental set-up, which can be 
divided into five major parts: 1) the RF section 2) mixer 
test dewar 3) dc measuring system 4) IF measuring system 
and 5) the computer-controlled instrumentation. In this 
chapter, we shall describe the experimental apparatus and 
discuss the calibration of the measuring system. 

5.2 The RF Section 


The millimeter-wave sources are two Varian VRT-2123A 
klystrons. One is used as the local oscillator, the other 
is used as a monochromatic signal source. The local 
oscillator frequency is set at 113-9 GHz while the signal 
klystron can be operated either at 115-5 GHz (upper 
sideband) or 112.5 GHz (lower sideband) so that the mixer 
performance at different sidebands can be measured 
separately. Two home-built 3 MHz-bandwidth phase-locked 
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loop systems, with an effective pull-in range of 200 MHz 
peak-to-peak [42], are used to stablize the klystrons. A 
thermocouple type power sensor (Anritsu model MP 82 B1 
with Ml 83A readout) and a frequency meter (Hitachi model 
2277) can be switched into the system to measure the 
output power of the klystrons and their frequencies. The 
modulator (Hughes model 45217 H-3000) following the signal 
klystron functions as a switch to control the signal power 
during the measurements. The signal and 10 power were 
combined in a high directivity 20 dB directional coupler 
(Baytron 3N-40/20) which is connected to the signal input 
waveguide of the helium cooled test dewar. The loss along 
the signal and 10 paths is carefully measured with less 
than 0.1 dB error (up to the RF input port of the test 
dewar). Since the Anritsu power sensor has a maximum VSVR 
OF 1.6 from 90 to 140 GHz and an accuracy of ± 3$ at 115-3 
GHz, obtained by calibrating against a calorimeter (TRG 
model 981/982), the overall uncertainity in the measured 
loss of this section is about 0.25 dB. 
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5 • 3 Cryogenic System - the Test Dewar 


Mixing experiments were performed in a modified 
Infrared Model HD-3 dewar, which has a 5" diameter liquid 
helium cold head surrounded by a 77 K radiation shield. 

The dewar, which holds 0.8 liter of liquid nitrogen and 
1.1 liters of liquid helium, can operate for more than 24 
hours when operated at 4*2 K and for about 8 hours when 
pumped to 2 K. The cold stage (in vacuum), including the 
mixer itself, a bias-tee and a home-made cold pad ( ^ 18 
dB), is usually maintained at 2 K by pumping on the liquid 
helium. The backshort settings can be adjusted from 
outside the dewar by a micrometer driven vacuum 
feedthrough unit which is coupled to the backshorts by two 
0.125" fiber glass rods. The mixer and dewar temperatures 
are monitored by two silicon diode temperature sensors 
(Lakeshore model DT-500 DRC). The intended function of 
the cold attenuator is four-fold: it reduces the room 
temperature blackbody radiation; it acts as a cold RF 
termination for the mixer with a very well defined 
temperature; and it provides the necessary mechanical 
support and heat sink for the mixer block. 

It is found that the attenuation of this cold 
attenuator, which is made by inserting a piece of wedge 
shaped kapton vane, metalized with nichrome, into the 
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copper waveguide block, varies slightly as a function of 
time, therefore it is necessary to measure its attenuation 
periodically. During the periods when these experiments 
were carried out, its attenuation, measured at liquid 
helium temperature, was found to be 17*8 ± 0.2 dB. The 
cold attenuator is connected to the room temperature 
microwave system via a 2.5" long stainless steel WR-10 
waveguide plated on the inside with approximately 
25 n-inch of copper, which is sealed by a room temperature 
WR-10 waveguide vacuum window. The window consists of a 
polyetheylene plug, epoxied into a short reduced-size 
waveguide section. The return loss of the window is 
greater than 20 dB over a 6 GHz bandwidth. The insertion 
loss of the entire input waveguide section, measured with 
the dewar cooled to liquid helium temperature, is 0.35 dB 
± 0.05 dB. The overall attenuation of the dewar input 
waveguide and cold attenuator is 18.15 ± 0.25 dB. 



5*4 dc Bias and Measuring System 


The SIS junction is dc biased through a liquid helium 
cooled bias-T which is also used to separate the IP ( 1 . 4 
GHz) signal from the dc bias. The advantage of the cooled 
bias-T is that it can properly heat sink all the leads to 
ensure that the junction reaches the temperature of the 
mixer block. The IP signal is isolated from the dc 
circuit by a 20 pP chip capacitor (ATC-175 series). A 
coil inductor, close to its parallel self-resonance at 1.4 
GHz, is used as an RF choke. The return loss of the 
bias-T is greater than 25 dB between 1 and 2 GHz. The 
mixer bias voltage is developed across a 50 ohm bias 
resistor (thin film rod type) in the bias-T by a current 
source external to the dewar. The dc I-V characteristic 
of the junction is measured by a conventional 4-lead 
methods: 2 wires monitoring the junction voltage and the 
other two measuring the voltage developed accross the 5 
ohm current sensing resistor. Two instrumentation 
amplifiers (INA-101CM) are used for current and voltage 
amplification. In order to protect the junction against 
static and transients, all bias and monitoring wires 
entering the dewar were filtered by 4-element low pass 
filters (SPRAGUE series 1JX2851). The insertion loss of 
the filter is 80 dB at 90 KHz with a low frequency roll 
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off of 25 dB/octave. The bias voltage 
at a preset value by a leveling loop, 
handy when tuning up the mixer. 


can be maintained 
This is extremely 



5 • 5 IF Measuring System 


The IF (1.4 GHz) output of the mixer, passing through 
the bias-T, was coupled to a room temperature IF 
radiometer/reflectometer system via two series-connected 
50 ohm stainless steel coaxial cables, which were jointed 
at the 77 K radiation shield to reduce the heat load to 
the cold stage. The IF test set, similar to that 
described by Weinreb and Kerr [78], was used to measure 
the IF output power of the mixer, its IF reflection 
coefficient and the attenuation of the dewar's IF section. 
The 1.4 GHz radiometer is calibrated against liquid 
nitrogen and room temperature matched loads so that the 
detected power can be read directly in absolute 
temperature (± 0.5 K). The mixer's IF mismatch is 
measured by injecting a 1.4 GHz noise signal through the 
directional coupler towards the output port of the mixer 
and measuring the power reflected. By comparing the 
reflected power (corrected for the round-trip loss in the 
the dewar's IF section) with the injected signal, the 
magnitude of the reflection coefficient at mixer's output 
port can be determined. 

In order to determine the performance parameters of 
the SIS mixer (L, T^ and IF-VSWR), it is necessary to know 
the attenuation of dewar's IF section and its effective 
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radiation temperature. This can be accomplished by two 
separate measurements as described below. 

Figure 5-2 illustrates the various noise 
contributions to the measured output noise temperature 
from the dewar's IF section, which is terminated by a cold 
50 ohm load at its mixer port. Assuming that the 
attenuation of the IF section is a TT , and its effective 
radiation temperature can be characterized by 
(towards the mixer port) and T £ (towards the room 
temperature end), the output noise temperature is given 
by: 


T. = T 


IF 


r I' 


T a1 11 


IF 


+ T 


a2 


+ (1 - 


I r I ) “ip 


(5-1 ) 


Where T r is the room temperature, is the bath 
temperature and r is the reflection coefficient between 
load and the IF cable. The output temperature, with a 
noise signal T^g injected into the IF cable, is: 


T 2 = ^ T r + T R/S } “IF* 1 ' + T a1 “IF + T a2 


+ ( 1 I r I ) a jp 


(5-2) 
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Replacing the cold load by short circuit and repeating the 
above measurements, we have (without the noise injection): 


T 3 " T r a IF + T a1 “IF + T a2 


(5-3) 


and with the noise injection: 


T 4 " ( T r + T n/s) “IF + T «1 “IF + T 


a2 


(5-4) 


From equations (5*3) and (5.4), the attenuation a^p is 
given by: 


“IP “ V (I 4 - V 7 T 


N/S 


(5-5) 


solving (5*1), (5.2) and (5*5) for |r| , we have: 


rr = (t 2 - T 1 ) / (t 4 - t 3 ) 


(5.6) 


knowing a^p and |r| , the effective radiation 
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temperatures, T^, and T q2 can be calculated from the 
following equations: 

T a2 = [T 2 * T 4 |r|2 - T B °IP 0 - l r l 2 )l 

/ (1 - |r | 2 ) (5.7) 

and 

T a1 = “ T r_ °IF " T a2^ ^ °IP 

With the mixer in place, similar measurements can he 
performed by biasing the junction at high impedance region 
and 50 ohm point respectively, also including the shot 
noise generated by the junction as an additional noise 
source (discussed in chapter 7)« The measured results are 
o Ip = 0.8 ± 0.05 dB, T q1 = 17 ± 5 K and T a2 = 20 ± 5 K. 
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5 • 6 Computer Instrumentation 


An Apple II plus (48 K) computer was integrated into 
the system to enable experimental data acquisition and 
analysis to be performed automatically. Whenever a good 
operating point was located (the tuning technique will be 
described in the next chapter), the mixer performance was 
measured as a function of junction's bias voltage, which 
can be controlled by the computer through a digitally 
programmable 12-bit D/A converter (Tecmar, Apple DAI 01). 

A home-built 2-channel Apple relay driver was used to 
control the RF signal source (via the modulator) and the 
IF noise during the measurements. The output of the 
radiometer is recorded by the computer via a 14-bit A/D 
converter (Tecmar, Apple AD211) and stored on diskette for 
use in the mixer performance analysis. Final results can 
be obtained from the printer. 
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CHAPTER 6. THE MIXING EXPERIMENTS AND RESULTS 


6.1 Introduction 


Functionally, a mixer is characterized by its signal 
and image conversion loss, equivalent input mixer noise 
temperature, and output impedance at the intermediate 
frequency. We have seen in chapter 2 that Tucker's 
quantum mixer theory can predict the performance of an SIS 
mixer. In this chapter we discuss how the performance of 
our SIS mixer, described in chpaters 3 and 4, is measured, 
and present the experimental results. We also compare the 
measured results with the predicted performance for our 
particular junction, using the 3-frequency quantum mixer 
theory. As shown in this chapter, the agreement between 
the experimental and theoretical results is very good. 

6 . 2 Experimental Procedure 


The SIS mixer block was mounted on a copper heatsink 
plate which is attached to dewar's cold head. The dewar 
case was first evacuated to a pressure of about 10 micron 
Hg and then precooled by liquid nitrogen to reduce the 
consumption of helium during the final cool-down. The 
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dewar was then cooled down to 4*2 K by transferring liquid 
helium. Usually the experiments are conducted at 2 K 
which is achieved by reducing the helium vapor pressure 
with a mechanical pump (Welch DUO-SEAL Model 1376). 

During the experiment, an external Helmholtz coil is used 
to generate a static magnetic field 50 gauss) to 
suppress the Josephson critical current to its first 
minimum. Optimizing the mixer performance at a given 
sideband, upper or lower, is quite straightforward. As 
pointed out by Rudner et al. [53], the mixer noise 
temperature at most of the operating points is very low, 
so it is only necessary to minimize the mixer conversion 
loss and avoid operating in unstable regions when tuning 
the mixer. Therefore, the dc bias voltage, input pump 
power and the backshort settings are all optimized to give 
the maximum readout of the IF radiometer when a signal is 
applied to the input. Since our IF frequency (1.4 GHz) is 
rather high and the mixer tuning is very frequency 
dependent, it is necessary to measure the mixer 
performance at both sidebands under the same operating 
conditions. The measurement method is as follows. 
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The IP mismatch is measured, as described in section 

* 

5*5, by injecting approximately 1000 K of noise towards 
the mixer IP port and measuring the power reflected. 
Correction is made for the loss of dewar's IP section. By 
measuring the IP temperature increment when a known RF 
signal power is applied to the mixer, the conversion loss 
can be calculated. The available conversion loss is 
obtained after correcting for the IP mismatch. The mixer 
noise can be calculated from the above measurements. 
Detailed calculations are given below. 


Pig. 6-1 shows the various contributions to the noise 

power measured by the 1.4 GHz IP radiometer at dewar's 

output port with A) no RF or IP noise power applied to the 

mixer B) a known RP power, T . or T. applied to 

sig image’ ■" 

dewar's RP input port and C) an IP noise power T N ^ g 

applied to dewar's IP port. Referring to this figure, the 

measured noise temperature, T tx1 , T T -n and T T -r, , in each 

IF A IF B IF C 

case is then given by: 


T IF a ~ ^ a RP T r + ^ “ a RP^ T B^ 1//L s + 1 / L i^ + T M //I s^ 


* All the noise powers mentioned in here are expressed in 
terms of the equivalent noise temperatures, T = P/kAv. 
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x [(1 - |r| 2 )a ip ] 


+ “l P m 2 [T a1 + a Ip T r ] + I a2 


( 6 . 1 ) 


IF 


= T 
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IF, 


+ T 


‘Big 


“rf ^ 1 
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“IP (1 - 


rl 2 ) 
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( 6 . 2 ) 
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T. 

image 


OrfO / li) a IF (1 


I r | 2 ) 


and 


T 




+ T 


N/S a IF 


I r r 


( 6 . 3 ) 


where 

a RF = a ^‘ ten ' ua ‘ti° n of dewar's RF section 
T r = room temperature 
= bath temperature 

r = reflection coefficient at mixer IF port 

L = conversion loss at signal frequency 
s 

L i = conversion loss at image frequency 
Tj^ = equivalent single sideband noise temperature of 
the mixer referred to the input port. 
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a IP* T a1 ’ and T a2 are define(i in section 5*5 as the 
dewar's IP attenuation and it's effective radiation 
temperatures, respectively. 

Solving (6.1) and (6.3) for |r| 2 , we have: 

I r I 2 = t 1 / »fpH(T Ip - T Ir ) / T n/s ] (6.4) 

C A 

From equations (6.1) and (6.2) the conversion loss at 
signal (image) frequency L s (L^ ) is given by: 


L s(i) " [ a RF a IF I ^ T IF “ T IF.^ 

£5 A 


o - in 2 ) 


sig( image) 


(6.5) 


Knowing , 1^ and |r| , the single sideband mixer noise 
temperature can be calculated from (6.1) as: 

- t T IF A - T „2 - “lFl r l 2 < T a1 + “ipV^S 

/ U IP (1 - I r | 2 ) ] 

- [a RF T r + (1 - « Rp )T B ][l + (L b / L.)] (6.6) 
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6.3 Experimental Results 


The mixing performances of various SIS arrays, 
fabricated by IBM and NBS-Boulder, have been tested at 115 
GHz, and are summarized in table 6-1. These results are 
discussed in detail below. 

6.3-1 The Type J-J Mixer 


Figures 6-2 and 6-3 show the typical experimental 

results for the type J-J mixer with the IBM 14-junction 

array. The dc I-V curve with or without 10 power (F^q) 

and the measured output noise temperatures T t1? , T x -p and 

1J? A 1X B 

T tp , which are described in the previous section, are 
plotted on the same dc bias voltage scale. Being operated 
at 2 K, the unpumped dc I-V curve shows a very sharp rise 
at the energy-gap voltage, V = 14(2i/e) a 41*5 mV, and a 

O 

barely noticeable subgap structure at about 23 mV. 

Measured at the center of the first photon step, V = V^ - 
14(hu/2e) = 38 mV, the leakage current is around 5 uA. 

The normal state resistance of array is rather high, R N = 
600 ohm measured at V = 55 mV. Since the array 
capacitance is 22.5 ± 2.0 fF, this yields an uR^C product 
of 10 for the entire array at 115 GHz. 




124 


With LO power applied, the pumped dc I-V curve 
exhibits the photon-assisted tunneling structure with a 
voltage step of Mliu/e. Similar voltage-modulated 
structures are also observed in the three measured output 
noise temperature curves and therefore in the mixer 
performance parameters (L c , T^ and IF-VSWR). The best 
conversion is obtained when junction is biased at the 
first peak below the gap voltage. This agrees with the 
general trends observed by Rudner et al. [53] that for 
large R^ the conversion is more efficient at the first 
peak while for low R^ the high number peak (low bias 
voltage) is more efficient. 

As shown in Fig. 6-4, under certain tuning conditions 

negative differential resistance is observed at the first 

photon peak. This implies that, for low IF, infinite 

conversion gain is available into an appropriate IF load. 

This has been verified by the dc incremental method [70]. 

The mixer's IF reflection coefficient, calculated from the 

T tp -V curve, is also found slightly greater than unity at 
C 

V = 38 mV. Therefore, the mixer's IF output resistance is 
negative and the infinite gain is also available at 1.4 
GHz. Since the array normal state resistance is extremely 
high, the mixer is very badly matched into the 50 ohm 
load. Therefore the maximum realized conversion in these 
experiments, from mixer input to the 50 ohm load, is about 
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11.5 dB. This is achieved "by biasing the junction at the 
center of the first photon step. By installing a high 
ratio IF transformer or by terminating mixer with an 
appropriate high resistance load, substantial conversion 
gain could be realized in this mixer. At the same bias 
point, the equivalent single-sideband mixer noise 
temperature is 70 ± 40 K. High noise due to the Josephson 
effect, which is sensitive to the applied magnetic field, 
was observed for bias voltages below 28 mV. At higher 
voltages, the mixer performance is independent of the 
applied magnetic field. The saturation power measured at 
the 1 -dB gain. compression point is about 30 nW. 

6.3*2 The Type-D Mixer 


Several NBS N-junction series arrays, (N =2, 4), 

were tested in the type-D mount. Some typical results are 

presented in Figures 6-5 and 6-6. The general features of 

the unpumped dc I-V characteristics of these arrays, such 

as the steep current raise at the gap voltage, V = 

§ 

N(2A/e), low leakage current and sharp corner below the 
gap voltage, are very similar to that of the IBM array. 

The normal state resistances of these arrays are between 
50 and 100 ohm. Using an estimated array capacitance of 
0.1 pF gives uR^C = 4 to 7* The pumped dc I-V curves show 
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very clear photon steps with voltage periodicity of Nhu/e. 
By adjusting the backshort settings and the input LO 
power, it is possible to obtain a negative differential 
resistance region on the pumped dc I-V curve at the first 
photon step. Measured under the same tuning conditions, 
an IF-VSWR about 1 .8 is obtained. The significant 
difference between dc and 1.4 GHz IF output impedance is 
due to the sharp RF tuning of the mixer. 

Since this mixer mount has very broad tuning range, 
the mixer can be tuned for operation with almost any 
desired sideband ratio. The performance of the 
two-junction arrays is generally superior to that of the 
4-junction ones. The best result obtained, with the mixer 
tuned for single-sideband operation (lower-sideband 
rejection > 25 dB) and biased at the first photon step, is 
a conversion loss of 6.9 ± 0.5 dB at 115-3 GHz and the 
mixer noise temperature of 20 K . The mixer’s 
instantaneous bandwidth is 350 MHz, and the 1 -dB gain 
compression point occurred at ^ 4 nW. 

This SIS mixer, whose performance is superior to any 
yet reported at 115 GHz, has been used to construct an 


* As shown in the next chapter, the mixer noise 
temperature deduced from receiver's performance is 15 ± 14 

K. 
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ultra-sensitive 110 - 115 GHz receiver. The construction 
and performance of this receiver will be described in 
chapter 7- 
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6 . 4 Comparison: Experimental Results Versus 
Theoretical Predictions 


Dr. M. J. Feldman has calculated the performance of 
our mixers by applying the three-frequency quantum mixer 
theory to the measured unpumped dc I-V curve of the array 
under the following assumptions: a) All of the junctions 
in the array are identical. Then the array can be treated 
as an "equivalent single junction", b) The junction 
capacitance is large enough to short circuit all of the 
higher harmonics, c) The dc current through the array is 
entirely due to the single particle tunneling. His 
results are used to compare with our experimental data. 

6.4.1 Quantitative Verification of Tucker's Theory f 1 8 ] 


We have carried out a series of measurements on the 
type J-J mixer (with the IBM 14-junction array) at 115 
GHz. The results are found to be in good agreement with 
the theoretical predictions. In these experiments, the 
mixer performance was measured, with the same applied LO 
power 180 nW) , at each of nine equally'-spaced backshort 
settings (every 5 mils). Figure 6-7 illustrates one of 
these results obtained with the backshort set to 0.248". 

In order to compare with the theoretical predictions, four 
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data points, as shown in Pig. 6-7, were chosen from each 
of the nine backshort settings at the first and the second 
photon steps. They are the conversion losses and the 
pumped dc currents at the bias voltage of 38 mV (= V 

& 

Nhu/2e) and 31 mV (= V - 3Nhu/2e). These 36 experimental 

O 

data points, displayed in Pig. 6-8 as a function of 
backshort settings, are used to compare with the 
theoretical calculations. 


Since the measured result of our mixer noise 
temperature, T^ = 70 ± 40 K, is not precise enough for 
quantitative comparison, no attempt was made to compare 
the mixer noise temperature. 

In order to simulate the experimental situation, 
Feldman used the following experimental information in his 
calculation: a) the experimentally measured incident LO 
power, P^q = 180 ±40 nW, b) the IP load resistance is 50 
ohm, c) the array capacitance, C = 22.5 ± 2 fP, and d) the 
complex embedding admittance at each backshort setting, 
for signal, pumped and image frequencies, measured on the 
scale model and shown in Pig. 4-17. 


One of the calculated results, with P-^q = 1 73 nW and 
C = 22.5 fP (both within the experimental errors), are 
shown in Pig. 6-8. The computed conversion losses and the 
pumped dc currents at 38 mV and 31 mV are displayed as a 
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function of backshort setting for comparison. The 
agreement between the theoretical results and the 
experimental data points is extremely good. This 
quantitatively verifies Tucker's three-frequency quantum 
mixer theory, at least for an SIS mixer with uR^C product 
about 10. This is the most convincing quantitative 
verification of Tucker's theory which has yet been 
published . 

It is found that good agreement between theory and 
experiment exists if the theoretical calculations assume 
P^q between 165 nW and 180 nW, and C between 22 fP and 23 
fF. The measured values of P^q and C are 180 ± 40 nW and 
22.5 ± 2.0 fF respectively. It is clear that, by fitting 
the predictions of quantum mixer theory to the 
experimental mixer measurements, the capacitance of an SIS 
junction can be determined more precisely than by other 
techniques [2,37,82]. 

Ve also compared our measurements with theoretical 
calculations which assumed equal terminations for signal 
and image frequencies. The observed huge discrepancies 
between the theoretical and the experimental results 
indicate that the embedding admittance at each RF 
frequency has a very significant effect on the mixer 
performance. Therefore, one should determine these 
quantities as accurately as possible. 
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In order to determine the effects of the nonlinear 
quantum reactance, the same data set was compared with 
calculations neglecting the nonlinear quantum reactance 
terms in eq. (2.58). The computed results still agree 
very well with the experimental ones. Thus, we conclude 
that the nonlinear quantum reactance has no detectable 
effect on our results. 

6.4*2 Analysis of the Type-D Mixer 


The theoretical performance of the type-D mixer, 
whose unpumped dc I-V curve shown in Pig. 6-5, was 
calculated by Feldman. For equal signal and image 
terminations, he predicted a best conversion loss of 1.5 
dB. But in our experiments, due to the sharp tuning of 
the mixer mount, only one sideband could be well 
terminated while the other was approximately 
short-circuited (see Figs. 4-14 to 4-16). Thus, we should 
compare our results with the predictions assuming shorted 
image case. Fig. 6-9 shows contours of the predicted 
conversion loss, assuming the image is short-circuited, as 
a function of the normalized signal source conductance and 
susceptance plotted on the Smith chart. These contours 
were calculated by assuming an IF load resistance of 50 
ohm and that the junction is biased at the first photon 
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step below the gap-voltage. Furthermore, at each of the 
signal source admittances, the junction bias voltage and 
the normalized pumped voltage a [= eV^/ (2ftu) ] are 
optimized for the best conversion. The optimum a is found 
between 1.2 and 1.6. As shown in Fig. 6-9, the best 
conversion loss, 4.3 dB, occurs at the optimum source 
admittance Y° pt = [(0.0165 ± 0.0002)-j (0.0045 ± 0.0002)] 
mho. Recalling that, in section 4*4.4, our simulations of 
the type-D mount embedding admittance have shown that in 
the vicinity of this optimum source admittance point, the 
circuit transmission loss will add 1-3 dB to the 
conversion loss. Hence in the short-circuited image case, 
the predicted overall conversion loss is 5-8 dB, which is 
in good agreement with our measured results. 

In the vicinity of the best operating point, the 
predicted mixer output noise temperature is 3-9 ± 0.4 K 
which corresponds to an equivalent input noise, temperature 
of 19 ± 4 K if the experimental conversion loss, 6.9 ± 0.5 
dB, is used in calculation. This agrees very well with 
the measured mixer noise temperature, T^ = 20 ± 15 K. 

As described in section 6.3«2, under certain tuning 
conditions, we observed a negative differential resistance 
region on the pumped dc I-V curve at the first photon step 
in this mixer. This phenomenon is also predicted by the 
theory. To show this, the predicted values of the local 


oscillator termination, which yields negative differential 
resistance on the pumped dc I-V curve, are mapped on the 
shaded area of the same Smith chart (Fig. 6-8). Since 
this shaded area falls within the available embedding 
admittance region shown in Fig. 4-9, we would expect 
regions of the negative dc differential resistance to 
appear on the pumped I-V curve. 
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Typical properties of lead alloy SIS array mixers operated at 2 K with f . = 115*3 GHz, 

f L0 = 115,9 GHz and f IP = 1,4 GHz 6 


TYPE OP 
ARRAY 

AREA OF EACH 
JUNCTION 
(nm 2 ) 

Jc 

(A/cm 2 ) 

r n 

(ohm) 

<jR n C 

L C 

(dB) 

I 1 1 
"nJ 

(dB) 

t m 

(K) 

IP- 

VSWR 

P L0 

(nW) 

SATURATION * ' 
POWER 
(nW) 

BAND-’ ' 

WIDTH 

(MHz) 

IBM 

14-JUNCTION 

8 

700 

600 

10 

1 1 .5±0.5 

-0.5 

70±40 

60 

160 

30 

- 

NBS 

2-JUNCTION 

8.4 

600 ° 

63 

8*5 

8. 0±0. 5 

>12 

20±1 5 

1 .05 

215 

7 

300 

NBS 

2-JUNCTION 

4*5 

760° 

94 

6.7 

6. 9±0. 5 

>25 

1 20±1 5 
'I5±14 n 

1.6 

250 

4 

350 

NBS 

4-JUNCTION 

17*2 

410° 

90 

12 

9.0±0.5 

- 

23±15 

1.5 

3400 

— 

— 


> l 
« i 
1 1 
« l 
1 1 


Image regection ratio, measured under the same tuning conditions 
Measured at 1 dB gain compression point 

1 dB bandwidth. This was limited by the mixer and not the IP system 
Assuming RmIq = 1.6 mV, typical for NBS junctions 
Deduced froni the receiver performance (see next chapter) 
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CHAPTER 7. A PRACTICAL SIS RECEIVER AT 1 1 5 GHz 


7-1 Introduction 


The previous chapters have explained the design, 
fabrication, measurement, and optimization of our SIS 
mixers, and have analyzed their performance. As shown in 
chapter 6, the agreement between the measured mixer 
performance and the theoretical predictions is excellent. 
It appears that the properties of our SIS mixers are well 
understood. The final part of this thesis is devoted to 
the development of a low-noise SIS receiver for 
radioastronomical applications at 115 GHz. In this 
chpater, we describe the construction of a 110 - 118 GHz 
ultra low-noise receiver using one of the type-D SIS array 
mixers which were optimized and described in chapter 6. 

We also discuss the performance of this receiver, which 
turns out to be the most sensitive one ever reported for 
this frequency range. 
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7-2 Receiver Description 


Pig. 7-1 is a simplified block diagram of the SIS 
receiver front-end. The local oscillator power is 
provided by a frequency-doubled 57 GHz klystron which is 
stablized by the same phase-locked loop described in 
chapter 5» The varactor diode doubler consists of a GaAs 
Schottky diode, mounted in a parallel WR-10/WR-15 
waveguide structure as described by Siegel [58], and can 
provide roughly 10 mV/ power at 114 GHz. The design of the 
receiver dewar is very similar to the test dewar described 
in chapter 5> the main difference being the use of a 
liquid helium cooled cross-guide coupler for 10 injection 
and the use of a liquid helium cooled field-effect 
transistor amplifier [79] as the first IP stage. Pig. 7-2 
is a picture of the complete SIS receiver. The dewar is 
an Infrared Laboratories' model HD-3(8) dewar which holds 
2.8 liters of liquid helium and 2-5 liters of liquid 
nitrogen and can provide approximately 1 1 hours continuous 
operation when the liquid helium is pumped down to 2.5 K. 
The incoming RP signal is collected by a room temperature 
scalar feed horn which is connected through a 3" long 
copper-plated stainless steel WR-10 waveguide and a 20 dB 
cross-guide LO injection coupler to the type-D SIS mixer. 
The total insertion loss from the aperture of the scalar 
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feed to the mixer is 0.6 dB and contributes about 12 K to 
the receiver noise temperature. The local oscillator 
power is introduced into the dewar via a 10" long 
stainless steel WR-10 waveguide to reduce the heat load to 
the liquid helium. Although the overall loss along the LO 
path is very high, around 30 dB, this scheme still 
provides enough 10 power to operate the SIS mixer. Both 
of the signal and the LO input waveguides are sealed by 
the room temperature waveguide vacuum windows described in 
section 5-3* 

The 1.2 - 1.6 GHz IP output from the SIS mixer is 
coupled through a bias-tee (described in section 5*4) and 
a 1 .4 GHz isolator (Passive Microwave Technology, model 
1102) to the L-band GaAs PET amplifier which is then 
connected to dewar* s room temperature IP output port via a 
0.085" stainless steel coaxial cable. As shown in Pig. 
7-2, all these major components are cooled to 2.5 K. The 
PET amplifier has about 30 dB gain and an input return 
loss > 1 5 dB from 1.2 to 1.7 GHz. The power dissipated 
from the amplifier is about 75 mW which is approximately 
half of the total heat load to the liquid helium. Thus 
any further reduction of amplifier's power dissipation 
would considerably increase the helium hold time of the 


dewar . 
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The insertion loss between the mixer and FET 

amplifier is 0.3 dB over the 1.2 to 1.6 GHz frequency- 

range which is mainly contributed by the isolator. The 

noise temperature of the IF amplifier can be measured by 

using the shot noise generated in the unpumped, dc biased 

* 

SIS array, as an IF noise source . Taking into account 
the mismatch between the mixer's IF output and the IF 
section, the noise temperature of the IF section is found 
to be 10.5 ± 1.0 K, referred to the mixer's output port. 
As described in section 5-5> all the bias and the 
monitoring wires entering the dewar were carefully 
isolated and filtered to avoid the RF pick-up and protect 
the junction against static discharges and transients. 

The mixer backshort adjustment assembly, temperature 
monitoring system and the junction bias and monitoring 
circuit are the same as those used in the test dewar. 
Details are given in chapter 5» 

The receiver was usually tuned for > 20 dB image 
rejection, thus operating as a true single sideband 
receiver. This was measured using a sideband generator 


^ As shown in [19]> if an N-junction array is biased with 
a dc current I. at a differential resistance R , , the shot 
noise power per c unit bandwidth is el^R^/^N. 
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which generates equal amounts of power at both the upper 
and the lower sidebands. The sideband generator consists 
of a 2.4 to 3*7 GHz varactor tuned transistor oscillator 
(Avantek, model VT0-8240) and a conventional GaAs Schottky 
diode mixer [12], which is used as a harmonic generator. 
The Schottky diode mixer was pumped at 2.812 GHz to 
generate 40th and 41st harmonics at 112.48 GHz and 115*29 
GHz respectively. The output at both sidebands was 
measured, the harmonic generator tuned so they were found 
to be within 1-dB difference. 

As shown in the previous chapter, the mixer 
performance depends on four strongly interacting 
parameters, namely: the dc bias voltage, the incident 10 
power, and the two backshort settings. Two independent 
leveling loops, as shown in Fig. 7-1, were used to 
simplify the mixer tuning process. The first leveling 
loop maintains the constant dc bias voltage which is 
usually set at the center of the first photon step below 
the gap voltage, the other one maintains the constant dc 
current through the junction (thus the constant pump 
voltage V L q across the junction through eq. (2.78)) by 
adjusting the incident LO power via a modulator (Hughes, 
model 4521 7H-3000) . The receiver performance can be 
optimized much more easily if these parameters are 
independently controlled. 
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7 • 3 Receiver Performance 


The receiver's response to room temperature and 
liquid nitrogen loads is shown in Fig. 7-3: the receiver's 
1 .4 GHz output power, in a 50 MHz bandwidth, and 
junction's dc current with and without the incident 10 
power 250 nV) are plotted on the same dc bias voltage 
scale. At the dc bias voltage of 5*45 mV, the measured 
single sideband receiver noise temperature, with an image 
rejection ratio > 25 dB, is 89 ± 3 K* at 1 1 5 • 3 GHz, 
referred to the aperture of the scalar feed horn. Since 
the junction was biased at the first photon step, no 
magnetic field was applied in these measurements. Fig. 

7-4 shows the various noise contributions to the overall 
receiver noise temperature at 115 GHz. Assuming a mixer 
conversion loss of 6.9 dB, as measured in the test dewar, 
the mixer noise temperature deduced from Fig. 7-4 is 15 ± 
14 K, which is in good agreement with the measured result 
presented in table 6-1 . 

The saturation power of the receiver is limited by 


T T-n has been improved to 83 K by installing a new SIS 
2-junction array and retuning the FET amplifier. The 
retuned amplifier dissipates less heat and has a lower 
noise temperature at the expense of a narrower bandwidth. 
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the SIS mixer whose 1-dB gain compression point occurs at 
4 nW input. As shown in Fig. 7-5> the receiver's noise 
temperature changes less than 10 K within a 270 MHz 
instantaneous bandwidth. 

If we operate the receiver at 4*2 K, the single 

* 

sideband receiver noise temperature increases to 200 K . 
This is mainly due to the higher mixer conversion loss, 
roughly 4 to 5 dB more, at 4.2 K. Properties of this SIS 
receiver are summarized in table 7-1 • 

Fig. 7-6 shows the single sideband receiver noise 
temperature as a function of signal frequency for an image 
rejection ratio larger than 20 dB, when it is retuned at 
each frequency. The SSB receiver noise temperature is 
about 80 K below 115 GHz and rises at higher frequencies. 
This increase in the noise is a result of the cross-guide 
10 injection coupler, which has an evanescent mode at 118 
GHz. Better results can be achieved by relaxing the image 
rejection, for instance, with 12 dB image rejection ratio, 
a single sideband receiver noise temperature of 68 K was 
obtained at 112.3 GHz. This receiver has been proved to 


7 At 4*2 K, Tp was improved to 1 43 K by installing a new 
SIS array and retuning the FET amplifier. 
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be very reliable and simple to use. It is presently in 
use on the Columbia-GISS 4-ft telescope. 
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Table 7-1 Properties of the SIS receiver at 2.5 K and 

4.2 K, with f = 115-3 GHz, f LQ = 113-9 GHz, 

and f ip = 1.4 GHz 


BATH TEMPERATURE 

2.5 K 

4.2 K 

SSB T r 

89 K [85 K] 1 ^ 

200 K [143 K ] 1 ^ 

hs 

6.9 dB 

— 

t m 

15 K 

— 

BANDWIDTH 

270 MHz 

— 

IMAGE REJECTION 
RATIO 

>• 25 dB 

> 25 dB 

TOTAL HEAT LOAD 
TO LIQUID HELIUM 

156 mV 

— 

LIQUID HELIUM 
HOLD TIME 

11 hr 

— 


1 ) Using the new SIS array and the retuned FET amplifier 



CHAPTER 8. 


OBSERVATIONS OP DISTANT MOLECULAR 


CLOUDS TOWARD THE CYGNUS-X REGION 
WITH A 4-FOOT MILLIMETER WAVE 
TELESCOPE 


8.1 Introduction 


The SIS receiver described in chapter 7 is now 
successfully in use on the Columbia-GISS 4-ft telescope. 
This receiver, which has a single sideband noise 
temperature of 83 K at 115 GHz, is currently the best in 
operation at this frequency and is some 25 times more 
sensitive than its room temperature predecessor. More 
than 20,000 CO spectra have been obtained in the 1983-84 
observation season with the new SIS receiver. The typical 
integration time per scan is only about 1 minute for an 
rms noise of 0.25 K with a bandwidth of 250 KHz per 
filterbank channel. A first order fit is generally 
sufficient for the baselines. 

In order to demonstrate the power of this new 
receiver, we, in collaboration with Y.-L. Huang, undertook 
a CO survey toward the Cygnus-X region. Results are 
presented in this chapter. 


8.2 The Columbia-GISS Sky Survey Telescope 


The Columbia-GISS 4-ft Sky Survey Telescope is 
located on the roof of the Pupin Physics laboratories at 
Columbia University in New York city. The telescope 
consists of four major parts: the antenna, the SIS 
receiver, the spectrometer, and the control system. These 
will be summarized in this section. 

The antenna is a Cassegrain reflector with a 
half-power beamwidth (FWHM) of 8.7 arc min and a beam 
efficiency of 0.75 at 115 GHz. It consists of a 4-foot 
diameter main reflector with f/D = 0.375 and a 7 inch 
secondary. The effective f/D ratio of the Cassegrain 
system is 3*79. The feed is a 5-inch long scalar feed 
horn with an aperture of 1.5 inches. For additional 
details of the antenna, see R. Cohen [8]. 

The front-end of the telescope is the liquid-helium 
cooled SIS receiver (usually operated at 2.5 K by reducing 
the liquid helium vapor pressure) described in chapter 7. 
The incoming radiation is collected by the scalar feed 
horn, downconverted to a 1.4 GHz intermediate frequency 
(IF) signal by the SIS mixer, and then amplified by the 
liquid helium cooled FET amplifier. The output from the 
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SIS receiver is further amplified by a low noise 
transistor amplifier with 50 dB gain and then down 
converted again to the second IF of 150 MHz. This second 
IF is amplified (by approximately 50 dB) and sent to the 
spectrometer. At 115 GHz, the complete system has a 
single sideband noise temperature, measured in front of 
the scalar feed, of 85 K with a 25 dB image rejection 
ratio. 


The spectrometer is a 256-channel filterbank similar 
to the standard NRAO design [59] • The spectral resolution 
is 250 KHz, which corresponds to a velocity resolution of 
0.65 km s -1 at 1 1 5 GHz, and a total bandwidth of 166 km s -1 . 

The output of each channel filter is detected by a 
square-law detector and then sent to an analog integrator 
which integrates the data for 48 msec. The outputs of the 
integrators are then passed to a 12-bit analog-to-digital 
converter for readout to the computer. 

An on-line minicomputer (Data General Nova 1200) is 
used for telescope guidance, data acquisition and 
preliminary data processing. Details of the telescope 
mount and drive system, the back end and the computer 
system are given in G. Chin [6] and T. Dame [ 1 5 ] • 



8 . 3 Observations 


All the spectra were obtained by position switching 
(i.e. by alternately observing the source position and two 
nearby CO emission-free comparison positions) at roughly 
15-second intervals, until the rms channel-to-channel 
noise was about 0.25 K. Typically, the integration time 
per scan was 1 minute. The receiver is calibrated against 
a room-temperature black-body by the chopper-wheel 
technique, which automatically compensates for changes in 
the atmospheric attenuation and expresses line intensities 
in terms of the antenna temperature. Correction for the 
beam efficiency yields the radiation temperature, T^, 
which is related to the radiation intensity, 1^, by the 

p 

Rayleigh-Jeans formula: I = 2kT^/A , where k is the 
Boltzmann constant and A is the wavelength. The filter 
bank was centered at -60 km s~^ with respect to the local 
Standard of Rest (from which the local galactic motion has 
been removed), and covered a velocity range from -143 km 
s" 1 to 23 km s" 1 . 

Our observations of the J = 1+0 rotational transition 
of CO at 115 GHz toward the Cygnus-X region were made in 
March 1984. We mapped a square region, roughly 4° on a 
side, centered at galactic coordinates sl • = 78.5°, b = 1° 
(Fig. 8-1). The map consists of 522 spectra and the total 
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integration time was about 13 hours. 

In this region, CO emission lines have been detected 
over a wide range of velocities, from -80 to 20 km s - ^ . 
Only the emission from molecular clouds which lie well 
beyond the nearby Cygnus Rift, will be discussed here. 
These clouds have radial velocities more negative than -20 
km s 

Seven well separated clouds have been identified 
whose velocities are between -80 and -20 km s - ^ . In Pig. 
8-2, we show a typical spectrum from each cloud. The 
spectra are remarkably flat so that it was sufficient to 
remove only a first order baseline from the raw data. 

Fig. 8-3 shows the spatial maps of the largest cloud. The 
CO emission is integrated over velocities from -55 to -35 
km s" 1 . 
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8.4 Discussion 


A previous search for CO emission in this direction 
was performed by H.-I. Cong [ll] using a 40-channel 
filterbank with 1 MHz resolution (corresponding to 
2.6 km s -1 at 115 GHz). Cong [ 1 1 ] detected no emission with 
an antenna temperature higher than 1 K between -50 and 
-25 km s . This was probably because spectral dilution 
smeared the emission. 

Approximate distances to the seven clouds discovered 
in the direction of the Cygnus-X region can be obtained by 
assuming that the orbits of molecular clouds about the 
galactic center ar circle. This technique is standard 
when stellar companions, whose distance could be more 
accurately determined, cannot be found for interstellar 
clouds. 

For galactocentric radii larger than that of the sun 
(R q ), a reasonable approximation is that the galaxy's 
circular rotation velocity is a constant 250 km s“^ . 
Accordingly, the radial velocity of a cloud which lies in 
the direction of galactic longitude i at galactocentric 


distance R is: 
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V r = 250 sinn [(R 0 /R) - l] 


( 8 . 1 ) 


where R q is the radius of the solar circle. By adopting 
the canonical value R q = 10 kpc and observing the cloud's 
radial velocities, we can solve for their galactocentric 
distances and then determine their distances from the sun. 
The kinematic distances of the seven clouds are estimated 
to be between 6.5 and 12.7 kpc, corresponding to 
galactocentric radii of 1 1 to 14.7 kpc. 

Cloud masses are estimated from the velocity- 
integrated CO intensity W(C0) (= J T R dV). Ve assume a 
mean molecular weight per molecule of 2.72 and an 
emperical conversion factor 2 x 10^ cm - K km“ s [55] 
between the molecular hydrogen column density and W(C0). 
The seven distant clouds are found to have masses between 
1 and 6 x 10^ M , typical of values found for giant 
molecular clouds in the inner Galaxy by Dame [13]* Table 
8-1 summarizes the physical characteristics of the seven 
distant clouds. 

Evidence for massive star formation, namely the radio 
HII region DR7 [48], is found only in the largest cloud 
(No. 5). Heavy foreground obscuration may have hampered 
the optical detection of young stellar objects near the 
clouds. Although embedded massive stars might be detected 
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by their radio recombination-line emission, no systematic 
search for such emission has been made in the Cygnus-X 
region. 

Giant molecular clouds are concentrated in the 
galaxy's spiral arms. Because so many large clouds lie 
within the small longitude range of our map, it is likely 
that our line of sight is tangent to a spiral arm. As 
shown in Pig. 8-4, this arm is possibly the extension of 
the well-known Orion Arm which may merge with the Perseus 
Arm and the Outer Arm in this direction. 

This survey clearly demonstrates the power of our SIS 
receiver. Compared with the original room temperature 
Schottky diode receiver, the overall system noise 
temperature (i.e., the noise contributed from both the 
receiver and the earth's atmosphere) has been improved by 
a factor of 5« Therefore, 25 times less integration time 
is required to produced a spectrum with the same signal to 
noise ratio. Furthermore, the quality of the spectra is 
excellent; only linear baseline removal is needed. Since 
atmospheric variability contributes significantly to 
baseline instability, the good baselines are at least 
partly attributable to the short integration time per scan 
required to obtain adequate signal to noise ratio. This 
performance clearly makes our receiver one of the most 
powerful tools currently available for use in 
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millimeter-wave astronomy and the ease with which 
optimization can he accomplished makes it very practical 
receiver. Many time consuming experiments, such as 
large-scale galactic surveys, will he made feasible or at 
least will benefit by this receiver's development. 



TABLE 8-1. PHYSICAL CHARACTERISTICS OF THE DISTANT MOLECULAR CLOUDS TOWARD THE 
CYGNUS-X REGION 


hk 

CENTER OF CLOUD 

PEAK 

err 

DISTANCE 
FROM SUN 
(kpc) 

07 

GALACTOCENTRIC 

DISTANCE 

(kpc) 

C2T 

EFFECTIVE 

RADIUS 

(pc) 

C JT 

MASS 

Oo\) 


< 4 ) 

El 

' VELOCITY 
(km a" 1 ) 

1 

79.8 

0.8 

-21 

6.5 

10.9 

66 

2.2 

2 

77.0 

0.0 

-21 

7.2 

10.9 

53 

1 .2 

3 

80.3 

1.5 

-32 

7.6 

11.5 

64 

1 .8 

4 

77.0 

0.0 

-32 

8.4 

11.5 

87 

2.9 

5 

79.0 

1 .0 

-42 

8.9 

12.1 

116 

5.8 

6 

78.5 

1.3 

-62 

11.1 

13.4 

83 

2.1 

7 

77.9 

0.8 

-75 

12.7 

14.4 

85 

3-9 


Notes: (1) Adopting a "flat" rotation curve with V = 250 km b " 1 and R Q = 10 kpc 

(2) (Area/O l/2 

(3) Adopting N(H 2 )/W(C0) = 2 * 10 20 cm" 2 3 K" 1 km " 1 b 
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CHAPTER 9. FUTURE WORK, SUMMARY AND CONCLUSION 


9 • 1 Directions for Future Work on SI S Receivers 


Although our SIS receiver is almost twice as 
sensitive as any previously reported, this result is still 
very far from the ultimate quantum-limited sensitivity 
predicted for an SIS receiver. The reasons for this are 
explained in chapters 4 and 6. Several experiments are 
being planned to overcome the present practical 
limitations and enable the SIS receiver to approach its 
theoretical performance. The main directions for future 
work are: 

1) All Superconducting Mixer Circuit. 

As indicated in Chapter 4, the backshort contact 
resistance and the ohmic circuit loss in the mixer mount, 
which is made of OFHC copper with a chromium-gold 
suspended stripline circuit, add 1 to 3 dB to the mixer 
conversion loss. It is hoped that we can reduce this loss 
by using superconducting circuits. Experiments using a 
niobium mixer block and niobium suspended stripline 
circuits are now under way. 
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2) Inductively Shunted Junction. 

With the help of a 500X scale model, we have designed 
a new SIS chip on which the junction is shunted by an 
inductive loop designed to tune out the large junction 
capacitance at the signal frequency. Our computer 
simulation indicates that this new design will provide 
better tunability: mixer conversion gain is expected along 
with a wider bandwidth ( 2 GHz). 

3 ) Niobium Tunnel Junction. 

Although lead alloy tunnel junctions have proved 
suitable for our present work, there are several reasons 
for considering refractory (e.g. niobium, tatalum) tunnel 
junctions. These junctions have higher electrical and 
mechanical strength, better thermal recyclability, and 
better long term stability. They also have higher 
critical temperature which should allow operation at 
temperatures achievable by present closed-cycle 
refrigerators, typically 4 - 6 K. 

Several niobium edge-junctions, fabricated by IBM 
Thomas J. Watson Research Center, have been tested. A 
double sideband conversion gain of 0.3 dB and a double 
sideband receiver noise temperature of 50 K have been 
observed in prilimary measurements. Although the junction 
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normal state resistance is low ( >/> 50 ohm in this case), 
this result is already comparable with that of the lead 
junctions. Increasing the junction normal resistance to ^ 
100 ohm is expected to give better results. 
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9«2 Ultimate Limits to SIS Receivers 


The success of this work indicates that in the near 
future millimeter-wave receivers using SIS mixers may- 
approach the ultimate sensitivity- limit imposed by quantum 
theory. The two outstanding questions still to be 
answered are: what high frequency limitations are there on 
SIS mixers, and what is the correct quantum noise limit 
for the SIS mixer? 

1) The High Frequency Limitations of SIS Mixers. 

Philips et al. [46] have reported SIS mixing results 
at 230 GHz. Recently, Sutton has observed SIS 
quasiparticle mixing, along with Josephson effect mixing 
for frequencies up to 388 GHz [65]. Based on the 
assumption that the high frequency limitation of an SIS 
quasiparticle mixer- is governed by the sudden onset of 
Josephson effect noise below the drop-back voltage 
[15,52,53,57], Tucker and Feldman have estimated the upper 
frequency limit for pure quasiparticle mixing, in lead 
alloy junctions, to be 250 GHz and 350 GHz for uR^C = 2 
and 10 respectively [76]. 

2) The Fundamental Limit to the Mixer Noise 
Temperature . 
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Theoretically, an SIS mixer should have a minimum 
mixer noise temperature of the order of fiu/k (= 5*5 K at 
115 GHz), the quantum noise [19»76]. The measuring 
technique described in this work was not precise enough 
(with uncertainty of ± 15 K = ± 3 "hu/k) to detect this 
noise. A new technique has been developed to measure T M 
with an accruacy better than ± 2 K, and it is hoped that 
we should be able to gain a better understanding of the 
practical noise limits of SIS mixers in the near future. 
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9*3 Summary and Conclusions 


The research reported in this thesis is in three main 
parts : 

1 ) An experimental SIS mixer was constructed and used 
to verify quantitatively the predictions of Tucker's 
3-frequency quantum mixer theory [18]. We also made the 
first observation of negative dc differential resistance 
and infinite available gain in an SIS mixer [33] > 
non-classical effects which are not seen in conventional 
diode mixers. 

2) Tucker's theory was then used to design an ultra-low 
noise SIS mixer for a 115*3 GHz spectral line receiver 
[ 43 , 44 ]. The noise temperature of this receiver, 83 K, is 
the lowest ever reported in this frequency range. The 
true single sideband nature of the SIS mixer removes the 
uncertainties usually associated with the unequal signal 
and image responses of conventional heterodyne receivers, 
and also reduces the overall system noise temperature by 
eliminating the sky noise at the image frequency. The 
receiver has proved to be very reliable and simple to use 
on the telescope. 
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3 ) The SIS receiver was used on the Columbia/GISS Sky 
Survey Telescope to detect seven distant molecular clouds 
towards the Cygnus-X region. 

In conclusion, we have achieved the original goals of 
this research project, namely, to study the performance of 
the SIS mixer at 115 GHz and to develop an ultra sensitive 
SIS receiver for use on the Columbia-GISS 4-ft telescope. 
Our results have demonstrated that the SIS mixer is likely 
to be the best candidate for sensitive millimeter-wave 
detection in the future. The success of this project may 
influence the choice of receivers for use in the Very Long 
Baseline Array (VLBA) and Millimeter Array, both planned 
by the National Radio Astronomy Observatory, the 
Submillimeter Array now under consideration by the 
Harvard/Smithsonian Center for Astrophysics, and other 
future millimeter wave radio astronomy facilities. 
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APPENDIX I THE EXPERIMENTAL DETERMINATION OF THE 
EQUIVALENT INDUCTANCES L^, 

AND L' 


The equivalent circuit of the inductance measurement 
setup (Fig. 4-6) is reproduced in Fig. AI-1 . 

The unknown inductance can he determined as follows. 
The power coupled to the 50 ohm load with finite shunt 
inductance L is: 

P L = [V 0 2 /(4 x 50)] x [4x 2 /(1 + 4x 2 )] (AI-1) 


where x = uL/50 is the normalized reactance and u is the 
signal frequency. Without the shunt inductance, the 
transmitted power is: 


P L=. = V 0 2 /< 4 * 50) 


(AI-2) 


Therefore, the ratio of these two powers is: 
= 4x 2 /(1 + 4x 2 ) 


r = V P L= 
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= 4(uI) 2 /[50 2 + 4(uL) 2 ] (AI-3) 

and the unknown inductance can be determined from the 
power ratio r by the following equation: 


L = (25/2»f)[l/(1/r - 1 ) 1/2 ] 


(AI-4) 


A 300X scale model (see Pig. AI-2) was constructed 
for inductance measurement. As shown in Pig. AI-3, all 
the inductance measurements were made relative to the 
inductance of the full width (W) stripline. The results 
are summarized in Table 4-1 • 



163 


APPENDIX II DETERMINATION OF THE EQUIVALENT CIRCUIT 
ELEMENTS OF THE TYPE-D MOUNT 


In this section, we will describe the experimental 
procedure used to determine the equivalent circuit 
elements of the two port network which represents the 
transistion between the junction and the single-mode 
waveguide. As illustrated in Fig. 4-11, this two port 
network can be represented by a T-network which is shown 
in Fig. AII-1 . 

We have used the same 40X scale model (described in 
section 4*3) to measure the elements of this T-network. 
The printed circuit pattern, as shown in Fig. AII-2, is 
modified so that it produces a short circuit at plane 
A-A'. The original filter pattern already produces a 
short circuit at plane B-B' at the signal frequency. 
Therefore in the equivalent circuit shown in Fig. AII-3, 
the series reactance Xg/g = 0 and the embedding impedance 
Z emb is ^ven by: 


Z emb = + IVIU& 0 + 1/(J‘ X l 2 + Z y/ G )]{ 

= f>Vc " V L 2 " V c + Z V/fi (dX l 1 + jx c )] 

/ (jX c + jX l2 + Zy/ G ) 


(AII-1 ) 
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It requires three independent measurements to determine 
the three unknowns: 

1 ) Measure the impedance looking into port 2 while 
shorting port 1 . Port 1 can he shorted by moving the main 
backshort to an appropriate position. The measured 
embedding impedance is: 

Z mea “ 3", = J X L, + 3 x L 2 X c ' < X I 2 + V < AI1 ' 2 ) 

2) Replace the broadband load in the input waveguide 
by a sliding short circuit and adjust for an open circuit 
at port 1. The impedance measured at port 2 is: 


Z mea = 3“2 = 3< X 1, + X o) (AII - 3) 

3) Adjust the impedance Zy/ G at port 1 until the 
impedance measured at port 2 is open circuited. This 
value of ZyjQ. is given by: 


Z W/G “ 3 m 3 " "^ X L 2 + X c) 


(AII-4) 
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Prom equation AII-2, 3> 4, we can solve for X T , X T and 

L 1 L 2 

X in terms of m. , m 0 and nu: 
c 12 3 


x c = -[m* x - m 2 )] 


1/2 


X L — m 2 + x — m2 ) ] ^ 


(AII-5) 


and 



■m3 + [m 5 x (m 1 



Measurements were performed in 25 MHz intervals from 
2.7 to 2.975 GHz corresponding to 108 to 119 GHz, the 
actual operating frequency range. Due to the residual 
errors within the network analyzer and the mechanical 
tolerances of the model, the accuracy of these 
measurements is limited to approximately ± 0.02 nH in L 
and ± 1 fP in C. The measured values of , L 2 and C are 
summarized in Table 4-2. 
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FIGURE CAPTIONS 


Fig. 2-1 Equivalent circuit of a mixer with large signal 
circuits to the right and small signal circuits to the 
left of the intrinsic junction which is characterized for 
small signals by the conversion matrix Y. The augmented 
network (broken line) includes all the sideband embedding 
admittances Y^ and is characterized for small signals by 
the augmented admittance matrix Y' . The equivalent signal 
current generator I is usually connected at port 1 of the 
augmented network during normal mixer operation, the other 
port being open circuited. 

Fig. 2-2 The mixer small signal equivalent circuit with 
shot noise current source and thermal noise current source 
at each sideband port. 

Fig. 3-1 Picture of a Princeton-NBS Pb(InAu)/oxide/Pb(Bi) 
2-junction array. Each junction is 1.8 um * 2.5 um, 
defined by the SiO window. 

Fig. 4-1 a) Type J-J mixer mount. It consists of an 
electroformed stepped waveguide transformer and a gold 
plated brass mixer body. A contacting backshort is used 
as tuning element, b) Side and top view of the SIS 



chip-Substrate assembly. The lead alloy 14-junction array 
was fabricated by IBM. c) Main body of the mixer mount. 
The substrate lies sideways in the slot and is held in 
position by the spring wire which also maintains the dc 
and RF/IF ground between the substrate and block. 

Fig. 4-2 a) Schematic diagram of the type-D mixer mount. 
The channel waveguide transformer is fabricated integral 
with the two halves of the block. The whole bolck is gold 
plated . 

Fig. 4-2 b) Top and side views of the SIS chip and quartz 
substrate. The SIS chip was mounted on the substrate by 
cold welding with indium. 

Fig. 4-3 Top view of the lower half of the 40X scale 
model of the type-D mixer block. Another coax cable, 
which is identical to the measurement cable, but short 
circuited at one end, is used as the phase reference for 
the network analyzer in the measurement. 

Fig. 4-4 A block diagram of the measurement setup used 
for the impedance measurement on the mixer model. 


Fig. 4-5 An enlarged view of the region in the vicinity 
of a) the SIS junction (quartz chip not shown) in the 
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actual mixer, t>) the corresponding region in the model 
(the upper block of Stycast not shown), and c) the 
reference planes A-A' and B-B'. The actual mixer and the 
model only differ within the reference planes A-A' and 
B-B'. Also illustrated are the equivalent circuit 
elements for the junction and the test probe inside the 
reference planes. 

Fig. 4-6 a) The measurement setup used to determine an 
unknown inductance, b) The equivalent circuit of the test 
setup. 

Fig. 4-7 The equivalent circuit used in deriving Y^^ 

seen by the SIS junction from Y which is measured on 

° mea 

the mixer model. 

Fig. 4-8 A Smith chart plot of the embedding admittance 
(Cj excluded) measured on the mixer model at 2.875 GHz 
(115 GHz in the actual mixer). Both conductance and 
susceptance are normalized to 1/(50 ohm). Each circle is 
obtained by fixing the side arm backshort at one position 
then moving the main backshort by half a guide-wavelength. 
The entire set of circles is obtained by repeating the 
whole process with different side arm backshort settings 
until it covers half a guide-wavelength. These circles 
cover all the admittance plane except for the forbidden 
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region indicated on the chart. The dotted line represents 
the corrected boundary of the forbidden region which 
includes the capacitive and inductive corrections shown in 
Fig. 4-5. 

Fig. 4-9 Available embedding admittance seen by the 
intrinstic junction at 115 GHz. A junction capacitance Cj 
= 99 fF ( uRjjC j = 6.7) is added in parallel to the external 
circuit . 

Fig. 4-10 An approximate microwave equivalent circuit of 
the type-D mixer mount. The circuit between the junction 
and the single-mode waveguide is characterized by the two 
port T-network. 

Fig. 4-11 A schematic drawing, showing the equivalent 
circuit elements surrounding the SIS array. They are 
idenfied as the elements of the T-network shown in Fig. 
4 - 10 . 


Fig. 4-12 A simplified microwave equivalent circuit of 
type-D mount. 

Fig. 4-15 Comparison of theoretically and experimentally 
determined forbidden region at 115 GHz. As described in 
the text, the measured values of and C, which are 
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small and have only */> 20$ accuracy, are used in the 
theoretical calculation. These errors caused the slight 
discrepancy shown in this plot. 

Fig. 4-14 One of the possible tuning conditions for the 
short-circuited image case, calculated from the equivalent 
circuit using the network analysis program described in 
section 4 . 4 . 4 . Here the embedding admittance (Cj 
included) seen by the intrinsic junction at 115*3 (+)> 
113*9 (■) and 112.5 GHz (•) is plotted on the Smith chart 
as a function of the two backshort settings. This plot 
was obtained by setting the main backshort at 0 . 058 " and 
then changing the side arm backshort setting, 0 . 0005 " per 
step, from 0.041" to 0.044". Also shown is the 
transmission loss (in dB) at each embedding admittance at 
the signal frequency. 

Fig. 4-15 One of the possible tuning conditions for the 
short-circuited image case, calculated from the equivalent 
circuit using the network analysis program described in 
section 4.4-4. Here the embedding admittance (Cj 
included) seen by the intrinsic junction at 115*3 (+)» 
113*9 (■) and 112.5 GHz (•) is plotted on the Smith chart 
as a function of the two backshort settings. This plot 
was obtained by setting the main backshort at 0 . 056 " and 
then changing the side arm backshort setting, 0 . 0005 " per 
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step, from 0.042" to 0.045". Also shown is the 
transmission loss (in dB) at each embedding admittance at 
the signal frequency. 

Fig. 4-16 One of the possible tuning conditions for the 
short-circuited image case, calculated from the equivalent 
circuit using the network analysis program described in 
section 4.4-4. Here the embedding admittance (Cj 
included) seen by the intrinsic junction at 115*3 ( + ), 
113*9 (■) and 112.5 GHz (•) is plotted on the Smith chart 
as a function of the two backshort settings. This plot 
was obtained by setting the main backshort at 0.0641" and 
then changing the side arm backshort setting, 0.0005" per 
step, from 0.034" to 0.037". Also shown is the 
transmission loss (in dB) at each embedding admittance at 
the signal frequency. 

Fig. 4-17 loci of the embedding admittance seen by the 
junction at the signal (S), local oscillator (10), and 
image (i) frequencies, measured on the 40X scale model of 
the type J-J mixer block. The array's capacitance Cj = 
22.5 fF has been added. 

Fig. 5-1 Microwave-cryogenic mixer testing set-up. 

Fig. 5-2 The various noise temperature components of 



172 


dewar's IF section terminated with a 50 ohm cold load. 

The whole section is represented by a attenuator with loss 

a IF ’ 

Fig. 6-1 The various contributions to the measured noise 
power at the dewar's IF output port with A) no RF or IF 
noise power applied to the mixer, B) a known RF power, at 
signal or image frequency, applied to dewar's RF input 
port, C) a known IF noise power, T^g, injected into 
dewar's IF port. The performance of the mixer can be 
determined from these measurements. a D1 , and a TT1 are the 
measured attenuations of dewar's RF and IF sections, 
described in chapter 5- 

Fig. 6-2 Typical experimental results of an IBM 
14-junction array mixer, a) The unpumped dc I-V curve and 
b) the dc I-V curve with a LO power of 160 nW at 11 3. 9 GHz 
applied to the mixer, c) - e) Shows the output of the 1.4 
GHz IF radiometer with c) a known IF noise power applied 
to the mixer's IF port, d) a 115*3 GHz monochromatic 
signal source (upper sideband) applied to mixer's RF port 
and e) no RF or IF power applied; all with P-^ = 160 nW. 

Fig. 6-3 Same as Fig. 6-2, except d) a 112.5 GHz 
monochromatic source was applied to mixer to measure the 
lower sideband performance. 
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Fig. 6-4 Similar experimental results for the same 

\ 

14-junction array mixer shown in Fig. 6-2 hut under 
different tuning conditions. The pumped dc I-V curve b) 
shows a region of negative differential resistance at the 
first photon step. The IF reflection coefficient, 
represented by curve c), is slightly larger than unity at 
38 mV. This implys a negative resistance at the 1.4 GHz 
IF. The structure seen in the curves below about 28 mV is 
very sensitive to the applied magnetic field and is caused 
by Josephson effect. The LO power is 375 nW. 

Fig. 6-3 Typical experimental results of an NBS 
2-junction array mixer, a) The unpumped dc I-V curve and 
b) the dc I-V curve with a LO power of 250 nW at 113*9GHz 
applied to the mixer, c) - e) Shows the output of the 1.4 
GHz IF radiometer with c) a known IF noise power applied 
to the mixer's IF port, d) a 115*3 GHz monochromatic 
signal source (upper sideband) applied to mixer's RF port 
and e) no RF or IF power applied; all with P^q = 250 nW. 

Fig. 6-6 Same as Fig. 6-5, except d) a 112.5 GHz 
monochromatic source was applied to mixer to measure the 
lower sideband performance. 


Fig. 6-7 a) The unpumped dc I-V curve for a 14-junction 
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SIS array, b) The dc I-V curve with P I0 " 180 nW applied 
to the mixer. c) The conversion gain with P^q = 180 nW 
and backshort position = 0.248". The four large dots, at 
38 mV and 31 mV, represent the experimental data points 
chosen from this backshort setting and used in Pig. 6-8 to 
compare with the theoretical calculation. 

Fig. 6-8 The four data points chosen from each of the 
backshort setting, as illustrated in Fig. 6-7, are plotted 
in here: a) the pumped dc current and b) the conversion 
gain, both measured at 38 mV with P^q = 180 ± 40 nW for 
nine backshort settings, c) and d) are the same as a) and 
b) respectively, but measured at .31 mV. The solid curves 
are the theoretical values calculated for C = 22.5 fF and 
P L0 = 173 nV. 

Fig. 6-9 Contours of the predicted conversion loss, 
calculated for optimum bias voltage (at the first photon 
step) and optimum a, assuming the image short-circuited, 
on the Smith chart of embedding admittance, normalized to 
1/(50 ohm). Negative differential resistance is expected 
to appear on the pumped dc I-V curve in the first photon 
step if the LO embedding admittance lies in the shaded 
region (again, a is optimized at each point). 

Fig. 7-1 Block diagram of the SIS receiver front-end, 
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showing the 10 source, receiver dewar, the LO leveling 
loop and the bias voltage leveling loop. 

Fig. 7-2 Picture of the completed 115 GHz SIS receiver. 
The principle elements on the 2.5 K cold stage are: 20-dB 
cross-guide coupler, SIS mixer, bias-tee, 1.4 GHz isolator 
and the L-Band FET amplifier. The 115*3 GHz RF input is a 
room temperature scalar feed horn. 

Fig. 7-3 a) dc I-V curve of the unpumped 2-junction SIS 
array at T = 2.5 K b) dc I-V curve with LO power of 250 n\J 
at 113*9 GHz applied to the mixer, c) and d) are the 
receiver's 1.4 GHz IF output with room temperature load ' 
and liquid nitrogen load placed in front of the horn. The 
receiver is tuned for > 25 dB image rejection ratio. 

Fig. 7-4 Breakdown of the SIS receiver noise temperature 
(SSB) at various points along the signal path (f . = 

O J. ^ 

115*3 GHz). T^ was deduced from the other measured and 
calculated quantities. 

Fig. 7-3 Reveiver's instantaneous bandwidth; data 
obtained with f^Q, P^q, c and backshort settings 
unchanged . 


Fig. 7-6 The receiver's noise temperature versus signal 
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frequency, for both upper-sideband (□) and lower sideband 
(0) operation. By installing a new SIS array and 
retuning the FET amplifier, T R = 83 K, indicated by (A) 
in this plot, was obtained at 115-3 GHz. The receiver is 
tuned for > 20 dB image rejection. Better results can be 
achieved by relaxing the image rejection ratio, e.g. with 
12 dB image rejection ratio, Tp = 68 K was obtained at 
112.3 GHz as noted on the figure. The uncertainty in T^ 
is ± 3 K. 

Fig. 8-1 Sampling map of the CO observation towards the 
Cygnus-X region; each observed point is marked by '+'. 

Fig. 8-2 Spectrum of the peak of each cloud listed in 
Table 8-1 . All spectra have a linear baseline removed. 

Fig. 8-3 Spatial map of the CO emission integrated over 
velocity from -55 "to -35 km s ^ . The contour interval is 
1.5 K km s" 1 . This is the largest cloud (No. 5) detected 
in this survey. 

Fig. 8-4 The locations in the galactic plane of the large 
clouds compiled in Table 8-1 . The circle sizes are 
proportional to the cube-roots of the cloud masses and a 
mass scale is shown. 
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Fig. AI-1 The equivalent circuit for inductance 
measurement. The unknown inductance is shunted across the 
50 ohm circuit and can be determined by measuring the 
power delivered to the load. 

Fig. AI-2 A 300X scale model used for inductance 
measurement . 

Fig. AI-5 a) The inductance (Lq) of the full width (V) 
stripline is used as the reference in the inductance 
measurements, b) and c) The magnetic energy stored in 
area I and II can be characterized by the equivalent 
inductances Lj and respectively, d) The inductance of 
the center conductor of the test probe. 

Fig. AII-1 The elements of the equivalent T-network 
between the junction and the single-moded waveguide. 

Fig. AII-2 A schematic drawing showing the modified 
printed circuit pattern used to determine the values of 
I 1 , I 2 and C. These circuit patterns produce a short 
circuit at planes A-A' and B-B'. 

Fig. AII-5 The equivalent circuit of the mixer mount with 
the modified stripline circuit pattern. 
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